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The macroscopic transport properties in a disordered potential, namely diffusion and weak/strong 
localization, closely depend on the microscopic and statistical properties of the disorder itself. This 
dependence is rich of counter-intuitive consequences. It can be particularly exploited in matter 
wave experiments, where the disordered potential can be tailored and controlled, and anisotropics 
are naturally present. In this work, we apply a perturbative microscopic transport theory and the 
self consistent theory of Anderson localization to study the transport properties of ultracold atoms 
in anisotropic 2D and 3D speckle potentials. We show that structured correlations can induce 
rich effects, such as anisotropic suppression of the white-noise limit and inversion of the transport 
anisotropy. We also calculate a disorder-induced shift of the energy states and propose a method 
to include it, which amounts to renormalize energies in the standard on-shell approximation. We 
show that the renormalization of energies strongly affects the prediction for the 3D localization 
threshold (mobility edge). We illustrate the theoretical findings with examples which are revelant 
for current matter wave experiments, where the disorder is created with a laser speckle. This paper 
is a detailed version of [Europhys. Lett. 99, 50003 (2012)] and represents a guideline for future 
experiments aiming at the precise location of the 3D mobility edge and study of anisotropy diffusion 
and localization effects in 2D and 3D. 
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I. INTRODUCTION 

Transport in disordered media is a fascinatingly rich 
field, which sparks a broad range of phenomena such 
as Brownian motion [1], electronic conductivity [2, 3], 
superconductivity [4], superfiuid flows of 4 He on Vycor 
substrates [5] , as well as localization of classical (electro- 
magnetic or sound) waves in dense media [6, 7] and of ul- 
tracold atoms in controlled disorder [8-12]. In the case of 
a matter particle for instance, two regimes should be dis- 
tinguished. In the classical regime, where the de Broglie 
wavelength is vanishingly small, transport leads to nor- 
mal or anomalous diffusion [13, 14]. The dynamics is 
characterized by the appearance of a percolation tran- 
sition, which separates a trapping regime - where the 
particle is bound in deep potential wells - from a dif- 
fusion regime - where the particle trajectory is spatially 
unbounded [15, 16]. In the quantum regime, the wave na- 
ture of the particle determines its transport properties, 
in close analogy with those of a classical wave [17, 18]. 
In this case, interference effects can survive disorder av- 
eraging, leading to striking effects such as weak localiza- 
tion [6], the related coherent back-scattering effect [19], 
and strong (Anderson) localization [20-22]. 

Localization is ubiquituous in wave physics. It shows a 
widely universal behaviour [23], but observable features 
significantly depend on the details of the system. It shows 
a renewed interest in the context of ultracold matter 
waves [8-12]. On the one hand, the microscopic param- 
eters in these systems are precisely known and, in many 
cases, tunable, which paves the way to unprecedented di- 
rect comparison between experiments and theory [24, 25]. 
This is a great advantage of ultracold atoms, compared 
to traditional condensed matter systems. On the other 
hand, these systems offer new situations, which can in- 
duce original effects [26] and provide new test-grounds 
in non-standard disorder [27-31]. Major advances in this 
field were the observation of one-dimensional (ID) An- 
derson localization of matterwaves [32, 33] and studies 
of the effects of weak [34-43] and strong [44-47] inter- 
actions in disordered gases. Presently, a major chal- 
lenge is the study of quantum transport in dimensions 
higher than one. While localization is the dominant ef- 
fect in one dimension [48, 49], higher dimensions show a 
richer phenomenology where regimes of diffusion, weak 
localization and Anderson localization can appear [23]. 
Recent experiments reported the observation of an An- 
derson transition in momentum space using cold-atom 
kick- rotor setups [50-52], study of classical diffusion in 
two-dimensional (2D) speckle potentials [53, 54], coher- 
ent back-scattering [55, 56], and evidence of Anderson 
localization in noninteracting Fermi [57] and Bose [58] 
gases in three-dimensional (3D) speckle potentials. 



From a theoretical viewpoint, diffusion and localiza- 
tion of noninteracting matter waves have been thoroughly 
studied for disordered potentials with zero-range correla- 
tions [59, 60] and isotropic correlation functions [61-66]. 

However, transport experiments in dimensions higher 
than one are most often done with speckle poten- 
tials which are anisotropic, either effectively in 2D se- 
tups [53, 54], or for fundamental optical constraints 
in 3D [57, 58]. Moreover, correlations in speckle po- 
tentials can be tailored in a broad range of configu- 
rations [67], which offers scope for investigation of lo- 
calization in nonstandard models of disorder [29, 30]. 
Taking into account anisotropic effects is of fundamen- 
tal importance because they can strongly affect coherent 
transport and localization properties. This was demon- 
strated for electrons in MOSFETs [68], diffusion- wave 
spectroscopy [69], biomedical imaging [70], and light in 
liquid crystals [71, 72], in phosphides [73], or in micro- 
cavities [74] . While the above systems are well described 
by models of disorder made with isotropic impurities em- 
bedded in anisotropic media [75, 76] or stretched scatter- 
ers embedded in isotropic media [77, 78], optical disorder, 
relevant to ultracold-atom experiments [57, 58], can show 
significantly more complex anisotropic correlation func- 
tions. The effect of the latter is much less known, and 
has been addressed only recently [79]. 

In this paper, we study quantum transport and Ander- 
son localization of matter waves in 2D and 3D anisotropic 
speckle potentials. We first introduce the basics of 
quantum transport of matter waves in disordered me- 
dia (Sec. II) and the models of disorder we focus on 
in 2D and 3D (Sec. III). We then present a detailed 
description of the theoretical framework pioneered in 
Refs. [75, 80], which intends to be pedagogical. Wc 
study single-scattering (Sec. IV), Boltzmann diffusion 
(Sec. V), and localization (Sec. VI), as a function of the 
particle energy, and discuss in particular the different 
anisotropics of these quantities. From a technical view- 
point, while the scattering allows for analytic expressions 
as for isotropic models of disorder [62], diffusion and lo- 
calization are more involved and require in general nu- 
merical diagonalization of a certain operator. Some an- 
alytic expressions are however found in some limits for 
anisotropic disorder. In Sees. IV, V and VI, we focus 
on the 2D case, which contains most of the anisotropy 
effects discussed in the paper. The 3D cases are dis- 
cussed in the next sections, where we study the same 
quantities as above (Sec. VII). We then show that energy- 
dependent quantities calculated in the usual on-shell ap- 
proximation should be renormalized in strong disorder, 
and propose a method to do it (Sec. VIII). Although it 
does not strongly alter the overall energy- dependence of 
the quantities calculated in the previous sections, it may 
be important when comparing to energy-resolved exper- 
imental measurements. Most importantly, we show that 
it strongly affects the calculation of the 3D mobility edge. 
Finally, we summarize our results and discuss their im- 
pact on recent and future experiments of ultra-cold atoms 
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Figure 1. (Color online) Schematic view of the coherent trans- 
port of a matter wave in a disordered medium, with special 
emphasis on the characteristic length scales. The figure shows 
a trajectory of a particle (solid multicolor line) in a two- 
dimensional disordered landscape (blue surface). Along its 
trajectory, the wave looses the memory of its phase (encoded 
in the various colors along the trajectory) on the characteris- 
tic length l s (scattering mean- free path). Multiple scattering 
then deflects the trajectory and the wave looses the memory of 
its direction on the characteristic length l B (transport mean- 
free path). Interference between the multiple-scattering paths 
can finally cancel diffusion (strong or Anderson localization). 
The wave then acquires an exponentially decaying probabil- 
ity profile (orange-green surface) of characteristic length L loc 
(localization length). 



in speckle potentials in the conclusion (Sec. IX). 



II. MATTER WAVES IN DISORDERED MEDIA 

A. Basics of quantum transport 

Before turning to a more formal description, it is worth 
recalling the basic picture of coherent transport in a dis- 
ordered medium. The basic ingredients are genuinely 
understood in a microscopic approach [80, 81]. Consider 
a wave of momentum k and velocity v propagating in 
a disordered medium [82]. Let us also assume, for the 
moment, that the medium is isotropic. We will drop this 
assumption in the following sections. The wave prop- 
agation is governed by scattering from the random im- 
purities. Three typical energy-dependent length scales 
can be identified [83], which characterize three basic ef- 
fects induced by the disorder (see Fig. 1). First, single 
scattering from impurities depletes the k-wave states, 
which can be seen as quasiparticles in the disordered 
medium, with a finite life-time r 3 (fc). Single scattering 
hence defines the first length scale, namely the scattering 
mean- free path, l a = vr a , which characterizes the typi- 
cal length travelled by the wave before it is looses the 



memory of its initial state, and primarily the memory of 
its initial phase. Then, multiple scattering defines the 
second length scale, namely the transport (Boltzmann) 
mean- free path, Z B , which characterizes the typical length 
travelled by the wave before it looses the memory of its 
initial direction. In general, several scattering events are 
necessary to significantly deflect the trajectories so that 
l B > l s . The two length scales are found to be equal only 
in the white-noise limit (if it exists), where the wave- 
length is smaller than the typical size of the impurities 
and the scattering is isotropic. In this case the wave 
looses the memory of its phase and initial propagation 
direction at the same time. Within the distance l B , the 
transport crosses over from ballistic to diffusive. The av- 
erage squared size of the wavepacket increases linearly 
in time, r 2 ~ 2D B t with D B = vl B /d the Boltzmann 
diffusion constant (d is the space dimension) [2, 3]. Fi- 
nally, diffusive transport allows the wave to return to its 
initial position via loop paths, and interference effects 
enter the game. Each loop can be traveled in one way or 
the other, which generates two multiple-scattering paths 
along which exactly the same phase is accumulated dur- 
ing the successive scattering events. This coherent effect 
holds for any specific realization of the disordered poten- 
tial and thus survives disorder averaging. Moreover, since 
these two paths are in phase, it gives rise to a construc- 
tive interference of the matter wave, which significantly 
enhances its return probability. This effect induces co- 
herent back-scattering and weak localization, which leads 
to diffusive transport with a reduced diffusion coefficient, 
£>* < D B [6]. For strong-enough disorder, the diffusion 
can completely cancel, an effect known as strong, or An- 
derson, localization [22]. Then, the probability distri- 
bution of the wave decays exponentially in space, hence 
defining the third characteristic length, L loc , the so-called 
localization length. 

The picture above shows that localization relies on 
two characteristics of the medium: coherence along the 
multiple-scattering paths and return probablity to the 
origin. One then understands that the strength of local- 
ization should be governed by the interference parame- 
ter kl B [84] (since the more the coherence length exceeds 
the typical length of a loop path, the more significant 
interference terms are), and by the dimension of space d 
(since the return probability decreases when d increases). 
As a matter of fact, in ID and 2D, any state is localized, 
although disorder correlations may lead to strong energy- 
dependence of the localization length [27, 28, 85, 86]. In 
ID, one finds that L loc ~ l B so that diffusion is strictly 
absent. In 2D, one finds l B < L lac , and diffusion shows up 
at intermediate distances and times. In 3D, the return 
probability is finite and localization appears only for suf- 
ficiently low values of kl B . A mobility edge shows up for 
kl B ~ 1, which separates localized states (for kl B < 1) 
from diffusive states (for kl B > 1) [23, 48, 87]. 

The microscopic description outlined above offers a 
comprehensive picture of transport and localization ef- 
fects for coherent waves in disordered media. In the 
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next subsections we will give mathematical support to 
this picture, and extend it using a formalism adapted to 
anisotropic disorder. 



B. Green functions 

Consider a quantum particle in a given homogeneous 
underlying medium and subjected to some static ran- 
domness. Its dynamics is governed by the single-particle 
Hamiltonian 



H = H„ 



V(t), 



(1) 



where Hq is the disorder-free, translation-invariant, 
Hamiltonian of the underlying medium, and V(r) is 
the time-independent (conservative) disordered poten- 
tial. For a particle in free space, which we will mainly 
consider below, the underlying medium is the vacuum, 
and H n = — ^;V 2 . Without loss of generality, the disor- 
dered potential can be assumed to be of zero mean value, 
V = [88] . The evolution of the wave function between 
to and t > to is determined by the retarded single-particle 
propagator 



G(t,t ) 



_ -iH(t-t )/h R 



e(t-to), 



(2) 



where the Hcaviside step function Q(t — to) accounts for 
temporal ordering. In the energy domain [89], G is the 
retarded Green operator 



G{E) 



[E -H + tiP 



(3) 



where E is the particle energy. It is the solution of the 
equation 

G(E) = Go (E) + Go (E) VG(E), (4) 

where Go = (E — Ho + iO + ) 1 is the disorder-free re- 
tarded Green function associated to the unperturbed 
Hamiltonian Ho- 

C. Properties of the disordered medium 

The calculation of any observable quantity is specific 
to the particular realization of the disorder. Therefore, 
meaningful quantities correspond to statistical averages 
over realizations of disordered potentials. When aver- 
aging over disorder realizations, some quantities can be 
written in terms of the average Green function G(E), 
for instance the spectral function (see below). The Born 
series of Eq. (4), averaged over the disorder, reads 



G = Go + GoVGo^Go + Go^Go^Go^Go 



(5) 



where the first order term Go^Go drops owing to our 
choice of energy reference, V = 0. It is convenient to 
represent this equation diagrammatically: 



(6) 



where a plain line is a Green function (grey for Go and 
black for G), the vertices (black dots) are scattering 
events and the dashed lines recall that they are corre- 
lated. One can write the Dyson equation [90] 



G — Go Go^jG, 



(7) 



where E(-E') is the self energy and can be developped in 
powers of V thanks to Eq. (5) [91]. The Dyson equation 
is formally solved by 



G = (Go 1 S) 



(8) 



Translation-invariant operators are conveniently writ- 
ten in momentum representation [92] , where they are di- 
agonal. For instance, the disorder-free retarded Green 
function reads 



(k|G o (£0|k') 



(27T) d J(k-k')Go(£,k) 

(2^(k-k ; ) 
E-e(k) + i0+' 



(9) 



where e(k) is the dispersion relation associated to Ho 
[(k|£Zo|k') = (27r) d (5(k- k')e(k)] and d the space dimen- 
sion. For an isotropic underlying medium, e(k) and the 
Green function Go(E, k) depend only on the modulus of 
k, k = |k|. If the disorder is homogeneous, i.e. if its 
statistical properties are translation- invariant [93], then 
the disorder-averaged Green function is also diagonal in 
k-space: 



(k|G(£)|k') 



(27r) d <5(k - k')G(E, k) 

{2n) d S{k - k') 
E - e(k) -£(£,k) +Ur 



(10) 



In addition, if the statistical properties of the disorder 
are isotropic, then G(E,k) = G(E,k). 

This features an effective homogeneous (i.e. 
translation-invariant) medium, which contains all 
necessary information to determine the disorder average 
of any quantity linear in G. It is the case of the spectral 
function A(E,k) defined by [81]: 



2Tr(k\S(E - H)\k') = (27r) d (5(k - k')A(E, k). (11) 

It contains all the information about the spectrum of the 
disordered medium. Using Eq. (3), it yields 



A(E,k) = -23 \G(E,k)] 



(12) 



The spectral function can be interpreted (up to a nu- 
merical factor) as the (normalized) probability density 
for an excitation of momentum k to have energy E and 
J ^-A{E, k) = 1. It is also the unnormalized proba- 
bility, per unit energy, to find a particle of energy E 
with momentum k and J A(E, k) = 2nN(E), where 
N(E) is the density of states per unit volume. For a 
particle in disorder- free space, it is given by A (E,k) = 
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k E =V2m E /h 



Figure 2. (Color online) Schematic representation of the spec- 
tral function A(E, k) of a particle of energy E, as a function 
of the particle momentum k. The vertical red line is the 
spectral function for the disorder free particle Ao(E, k) = 
2ir6 [E — e(k)] with e(k) = fi 2 fc 2 /2m. In the presence of dis- 
order the spectral function is shifted and broadened (black 
line). The standard on-shell approximation consists in ne- 
glecting the energy shift related to the real part of the particle 
self energy and the structure of the spectral function (dashed 
blue line). 



2ttS [E — e(k)]. In the presence of disorder, Eqs. (10) and 
(12) yield 



A(E,k) 



-2E"(£,k) 



(£-e(k) -£'(£, k)) +£"(£, k) 2 ' 



(13) 



with £' and S" the real and imaginary parts of S re- 
spectively. As represented schematically in Fig. 2, for a 
particle in free space [e(k) = h 2 k 2 /2m] with a weak disor- 
dered potential k) weakly depends on the momen- 
tum], the spectral function has a Lorentzian-like shape 
as a function of k. It is centered in ko, solution of 
E - e(ko) - £'(£,ko) = 0. The quantity £'(£,k ) thus 
describes the shift in energy of the free-particle modes 
when they are dressed by the disorder. The quantity 
H"(E, k) is the energy width of the spectral function, 
which defines the scattering mean free time 



r*(E,k) 



2T,"(E,k)' 



(14) 



or equivalently the scattering mean free path l s (E,k) = 
\v\T s (E,k). It accounts for the depletion of the free par- 
ticle mode at E = e(k) due to scattering from the disor- 
dered medium. Therefore the spectral function contains 
all the information about the relative weight, the energy, 
and the lifetime of the quasi-particles, i.e. the particles 
dressed by the disordered medium, which on average de- 
fine an effective medium. 



The spectral function will be a key element at each 
step of the following calculations, where we will compute 
quantities that depend on the energy E. In addition, in 
ultracold atomic systems, a broad range of energies are 
involved, but only the momentum distribution is usu- 
ally measured by time-of-flight techniques. The spectral 
function relates the energy distribution (T> E ) and the mo- 
mentum distribution (T> k ) of the stationary particles in 
the disorder via 



f dk 

V*{E) = J A(E, k)V k (k), 



(15) 



which is normalized by J ^-T> E (E) = 1. The exact cal- 
culation of the spectral function requires the knowledge 
of the real and imaginary parts of the self energy S [sec 
Eq. (13)], or, according to Eq. (11), the direct diago- 
nalization of the disordered Hamiltonian and an average 
over disorder realizations. This is, in general, a compli- 
cated task, especially in dimensions larger than one and 
for anisotropic disorder. In the following (see part IV) 
we will calculate the self energy by employing a pertur- 
bation theory and retaining only the first order (Born 
approximation) in Eq. (5) or (6). In Sees. V to VII, 
we work within the usual on-shell approximation [80], 
in which one neglects the real-part of the self energy 
J^'(E,k) and the structure of the spectral function (see 
schematic dashed blue line in Fig. 2). In Sec. VIII, we 
describe a method to go beyond the on-shell approxima- 
tion, which amounts to renormalizing the energies in a 
self-consistent way [79]. 



D. Propagation of the Wigner function 

Some quantities are not simply related to the averaged 
Green function G and require a more elaborate treat- 
ment. It is for instance the case of the spatial density 
and the momentum distribution. More generally, con- 
sider the time evolution of the one-body density matrix 
p(t) [81] or equivalently of the Wigner function [94] 



Pit) 



k - 



(16) 



The spatial density probability is given by n(r, t) = 
I (2Tr) d ^( r i t) and the momentum distribution by 
£> k (k, t) = JdrW(r,k, t). It is fruitful to rewrite 
Eq. (16) in a form indicating explicitly the initial con- 
ditions, making use of the relation p(t) = <d(t — 



t )e 



-iH(t-t )/fr 



p(t )e 



+iH(t-t )/h 



When averaging over 



disorder, if there is no correlations between the initial 
state and the disorder, one finds [63] 

W(r,k,t) = J dr' J ^W (r',k')F k y(r-r';t-t ), 

(17) 

where Wo(r, k) = W(r, k, to) is the initial Wigner func- 
tion and i^k.k' (R; t) is the phase-space propagation ker- 
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nel, defined by (if t > 0) 



and 



dE 



dq 



dfku 
~2^~ 



(18) 



<k + |G(£ + )|kV)<k'_|Gt(£_)|k_ 



(19) 

= (27r) d <5(q-q')<i> k , k '(q^,ii;), 

with k± = k ± q/2, k' ± = k' ± q'/2, E ± = E ±frw/2, 
and (q, w) the Fourier conjugates of the space and time 
variables [95] . As discussed above, disorder averaging fea- 
tures a translational invariance in space and introduces 
an effective medium for the expanding wave. As a di- 
rect consequence. Eq. (18) depends only on the differ- 
ence R = r r', which expresses the equivalence of all 
points in space. For the same reason, translational invari- 
ance, or equivalently momentum conservation, imposes 
that the sum of the in-going wavevectors (k + and k'_) 
on one hand, and out-going wavevectors (k', and k_ ) on 
the other hand, are equal. It leads to the condition on 
momentum transfer: q = q' in Eq. (19). 

As can be seen in Eqs. (17) and (18), the building 
block to describe wave propagation in random media is 
the density propagator $, which can be represented as 
a four-point vertex with k± and k' ± the left and right 
entries [see left-hand side of Eq. (21)]. The skeleton of 
this vertex is made by a retarded and an advanced Green 
functions (respectively G, represented by the top line, 
and G^ , represented by the bottom line). It contains 
all possible correlations between the scattering events of 
these Green functions. Following the same approach as 
used for the average field propagator G [leading to the 
Dyson equation (7)], the vertex <3? = G (g> Gt is formally 
constructed from the uncorrelated-average vertex G<E)G^ . 
Without any approximation, $ is then governed by the 
so-called Bethe-Salpeter equation (BSE) [81] 



$ = G®Gt + G®GtU$ 



represented diagrammatically as 
k+ k' + k+ 



k'_ 









+ 


u 











k' 



(20) 



(21) 



k' 



where U is the vertex function including all irreducible 
four-point scattering diagrams: 



U 




(22) 



The first term in the BSE [Eq. (20) or (21)] describes 
uncorrelated propagation of the field and its conjugate 



in the effective medium. The second term accounts for 
all correlations in the density propagation. 

Analogously to Eq. (8), the solution of the BSE (20)- 
(21) can be formally obtained from the inverse, if it ex- 
ists [96], of the four-point operator A = l-G(g>Gt U [97]: 



$ = A" X G® Gt. 



(23) 



More explicitly, the (k, k') component of a four- 
point vertex A which fulfills momentum conserva- 
tion is A k k / (q> w : E), such that (k + , k'_|A|k' + , k_) = 
(27r) d <5(q-q')A k , k '(q,a;,£), and 

A k , k ,(q,w,£0 = (27r) d «5(k-k') 

-Mq.w.^C/i^q.w.tf), (24) 

and 

/ k (q, u, E) = G(E+, k+)Gt (£L, k_). (25) 
Therefore Eq. (23) reads 

*k,tf(q,w,£) =A- 1 k ,(q,o; ) £;)/ k /(q ) u; ) .E), (26) 
and can be expressed as a geometric series 

$ k , k ' (q, w, E) = (2ir) d S(k - k')/ k (q, u, E) 

+/ k (q, u, E)U K u (q, lu, E)fa (q, u, E) 

+ J ^y I / k (q,a;, J B)C/ k>kl (q,a;, J E;)/ kl (q,a; ) S) 

xUk u u'(<l,u,E)fv(ci,Lj,E) + ... 

Note that the operator A _1 (cl),£') can be expressed in 
terms of the eigenvectors and associated eigenvalues of 
the operator A(u>,E). The calculation of the eigenvec- 
tors of the operator A(w, E) is the basic idea followed in 
Refs. [75, 98] to solve the BSE. It then gives access, via 
$, which is the quantity of interest, to the time depen- 
dence of the Wigner function [see Eqs. (16) to (18)], and 
of the spatial density in particular. 

In the following we will see that the intensity kernel $ 
has a diffusion pole, which takes the form 



$ k , k <(q,w,£;) 



A{E,k)A(E,k') 



2ttN(E) ifrui- nq-D(w,J5)-q 



(27) 



where D is the so-called dynamic diffusion tensor. The 
average spatial density distribution is then given by 



n(r,t) 



— W(r,k,t) 



(27T) 



= /S / dr ' P o(r',£)P(r-r',t-< |£l28) 

where V (r',E) = J -0^ A(E, k')W (r', k') represents 
the initial joint position-energy density and P(r — r', t — 
to\E) is the probability of quantum transport, i.e. the 
probability for a particle of energy E originating from 
point r' at time to to be in r at t. It can be expressed 
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thanks to Eqs. (17), (18) and (27) as the space-time 
Fourier Transform of the diffusion pole 1 / [ihu — hq ■ 

In the long-time limit, we will encounter two differ- 
ent situations. First, if \im ul ^o'D(uj, E) = T)(E) is a 
real definite positive tensor, the diffusion pole of the in- 
tensity kernel (27) describes normal diffusion with the 
anisotropic diffusion tensor D(E). We then find 



P(R,t -> oo\E) 



e -R-D _1 (_E)-R/4t 

^(47rt) d det{D(£0} 



0(t). (29) 



Second, if B(ui,E) - 0+ - iuA(E) in the limit uj ->• 0+ 
with A(E) a real positive definite tensor, the pole de- 
scribes localization. It leads to exponentially localized 
phase-space propagation kernel and probability of quan- 
tum transport at long distance. In 2D, 



P(R,i-> oo\E) 



e(t) (30) 



2^det{L loc (i;)} 
where Kq is the modified Bcsscl function, and in 3D 



P(R,i-> oo | £) 



^-v/R-L-^BJ-R 



4^det{L loc (£;)},/R • L^(E) • R 



=e(t). 

(31) 

In both 2 and 3D, the fonction -P(R) decays exponen- 
tially [99] over the characteristic length c (E) along the 
eigenaxes u of the localization tensor 1i\ oc (E) = yj A(E). 



E. Conductivity and Einstein's relation 

Finally, another quantity of interest for our problem 

- in parallel of those studied in sections II C and II D 

- is the conductivity. In complete analogy to the usual 
conductivity of charge in condensed matter systems [90] , 
we here define the conductivity tensor <x in our system as 
proportional to the current-current correlation function, 
via the Kubo formula [6, 100]: 



i (u,E) = 



(32) 



dk dk' 



/ (Zy (gp * [^<k|G(£? + )|k>^|Gt(g- 

where vi — hki/m is the velocity along axis i. As the 
structure of Eq. (32) is reminiscent of the definition of 
the four-point vertex $ [see Eq. (19)], calculations of the 
conductivity tensor can also be represented diagrammat- 
ically. The skeleton diagram, shown in Eq. (33), consists 
of the in and out-going velocities v and v' and of a bubble 
made of a retarded (top line) and an advanced (bottom 
line) Green function. As for $, the scattering events of 
the top and bottom lines can be correlated [see for ex- 
ample Eqs. (21) and (22)]. 



(33) 



Thanks to Einstein's classical argument, it was real- 
ized that, at thermal equilibrium, in a gas submitted to 
a force, the diffusion and drift currents have to be equal. 
This relation holds in general for quantum systems in the 
linear response regime (see e.g. Ref. [81]). In particular, 
here we expect the DC conductivity and diffusion ten- 
sors to be proportional : <t{lj = 0) oc D. More precisely, 
calculating cr B (u; = 0) in the Boltzmann and Born ap- 
proximations for anisotropic disorder permits us to find 
the proportionality factor (see details in appendix CI). 
In our system, we have 



(34) 



III. DISORDER CORRELATION FUNCTION 

Having recalled the general theory of quantum trans- 
port in disordered media, we now specify the frame- 
work of our study. In the following, we will consider 
ultracold matter waves as realized in several experiments 
[32, 46, 47, 53, 54, 57, 58, 101-106]. In our case, the 
underlying (disorder-free) medium is the vacuum, and 
e(k) = h 2 k 2 /2m. 

The models of disorder we will consider belong to the 
class of speckle potentials, which are particularly suited 
for ultracold atoms for they can be controlled [10, 12, 67]. 
In brief, a speckle pattern is created when a coherent 
light beam is shone through a diffusive plate (D) and fo- 
cused by an optical lens of focal distance / (see Fig. 3 
and Ref. [107]). At each point of its surface, the diffu- 
sive plate imprints a random phase on the electric field. 
The resulting electric field in the right-hand side of the 
lens is then the summation of many complex indepen- 
dent random components, and is therefore a Gaussian 
random variable according to the central limit theorem. 
The potential acting on the atoms is proportional to the 
intensity pattern (i.e. the square modulus of the electric 
field). It is thus a spatially random variable, but it is not 
Gaussian. However, thanks to the underlying Gaussian 
process for the electric field, the correlations of V(r) at 
any order can be obtained from the two-point correla- 
tion function of the electric field [28]. Here we consider 
weak disorder and, as we will sec, the two-point correla- 
tion function of V(r), C(r) = V(r)V(0), is sufficient to 
characterize the disorder in this case (we recall that we 
have chosen the zero of energy such that V = 0). 

For a fine-grain diffuscr, the two-point correlation func- 
tion C(r) is determined by the pupil function J D (p) (i.e. 
the intensity pattern just after the diffusive plate) [107]. 
In the following we will consider anisotropic Gaussian 
laser beams of waists w XtV and plates with homogeneous 

transmission, so that I D {p x ,p y ) = I oe - 2 ^ w ' + pV w l\ 
For the configuration of Fig. 3, in the paraxial approxi- 
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• ground 
glass 
plate (D) 



and we recover the power spectrum of a white noise dis- 
order, the only relevant parameter being V^a± x <T± y . 

The power spectrum (37) is obtained by shining an 
anisotropic Gaussian beam on the diffusive plate. It also 
approximately holds in the case of Ref . [53] where a quasi- 
2D Bose gas of width l z is subjected to a speckle created 
by an isotropic Gaussian laser beam shone with an angle 
9 with respect to the plane of atoms, if l z <C cr± <C ctii . 
In this case £ ~ 1/ sin6> (9 ~ ir/6 for the experiment of 
Ref. [53]). 



Figure 3. (Color online) Schematic of the apparatus used to 
create an optical speckle pattern. A laser beam is diffracted 
by a ground-glass plate diffuser (D) of pupil function I-d(p), 
where p = (px,p y ) spans the diffuser, which imprints a ran- 
dom phase on the various light paths. The intensity field, 
I(r), observed in the focal plane of a converging lens, is a 
speckle pattern, which creates a disordered potential V(r) for 
the atoms. 



mation, we find 



with 



ci sp (x,y,z) 



C(r) = V 2 c lsp (x,y,z) 



(35) 



exp 



2 / 2 



l+4z 2 /<7|| 



exp 



y 2 /<?± 

' l+Az 2 jo 



V 1 + 4 * 2 K 

(36) 
where A L 



TTW 



x.y 



a \\x, y = /^ w i,y and <r±x,v = 

is the laser wavelength. Here x and y are the coordi- 
nates orthogonal to the propagation axis z, and z = 
corresponds to the focal plane. It should be noted from 
Eq. (36) that the correlation properties in the (x, y) plane 
strongly depend on the value of z. We have chosen 
V n = y/C(r = 0) as definition of the amplitude of the 
disorder. 

In the following we will consider various cases in 2D 
and 3D that all derive from this configuration. 



B. Single speckle (3D) 

In Ref. [57] a 3D disorder is obtained by a single 
isotropic Gaussian laser beam of waist w (so-called single- 
speckle configuration). The disorder correlation function 
C(r) is given by Eq. (35) with w x = w y = w. The re- 
sulting speckle pattern is significantly anisotropic. It has 
correlation lengths a\\ in the propagation axis (z) and o~± 
in the orthogonal plane (x,y). In general 4/ > w, and 
C(r) is elongated along z (for instance <t\\/<j± — 5.8 in 
Ref. [57]). The corresponding disorder power spectrum 
reads [951 



with 



ci S p(k) 



3/2 _ 



C*(k) 



Va 2 Clsp(k) 



(38) 



Z-2 _i_ U2 



(39) 

It is isotropic in the (k x , k v ) plane but has a significantly 
different shape along the k z axis. This can be seen in 
Fig. 4(a): A typical iso- value surface of C(k) is wheel- 
shaped, and its cuts in planes containing h z are shaped 
like "8" figures. It also shows a strong algebraic diver- 
gence when k z = and k x + ky — > 0. It features absence 
of white-noise limit, which reflects the long-range correla- 
tions of the potential. The consequences of this property, 
obtained in the paraxial approximation, will be further 
discussed in the following. 



A. Anisotropic Gaussian speckle (2D) 



C. Orthogonally crossed speckles (3D) 



If the atoms are confined in a 2D geometry by a strong 
trapping potential along z centered on z — 0, they expe- 
rience a disordered potential with correlation function 



C(x,y) 
with a 



V£c lsp (x,y,0) = V^exp 



l Ax 2 + ey 2 ) , 

= a± x and £ = a± x /a± y the configuration 
anisotropy factor [108]. The Fourier transform [95] gives 
the power spectrum 



C(k) = V^V^exp 



V 



+ ^2") 



(37) 



Without loss of generality, we assume in the paper that 
f > 1. When |k| < OL-Vxy 1 , we get C(k) ~ V*ir?£ 



In Ref. [58], the disorder results from the crossing of 
two orthogonal speckle fields, propagating along the z 
and x axes, respectively. The polarizations and the fre- 
quencies of the two lasers can be chosen such that the two 
beams are fully incoherent (so-called incoherent- speckles 
configuration) or mutually coherent (so-called coherent- 
speckles configuration, case of Ref. [58]). 



1. Fully incoherent case 

In the incoherent-speckles case the 3D-correlation 
function in real space and the corresponding power spec- 
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Figure 4. (Color online) Disorder power spectrum C(k) for the (a) single-speckle, (b) incoherent-speckles, and (c) coherent- 
speckles cases [Eqs. (38), (40), and (41)] with the parameters of Refs. [57, 58] (<j\\/<t±_ — 5.8, and for (c) A L /o"x — 2.16). 
The functions C*(k) are represented as iso- value surfaces (at 2V^a\) and cuts in the planes defined by the symmetry axes: 
{u x ,u y ,ii z } for (a) and {u x = (u^-UzJ/V^.uy = u y ,uz = (u x +u z )/y/2} for (b) and (c). 



trum C(k) are given by the sum of two orthogonally- 
oriented spectra, similar to that of the single-speckle case, 
so that 

C(k) = (14/2) 2 [ci ep (k x ,k y) k z ) + ci sp (k Zl k yi k x )] , 

(40) 

where 5i sp (k) is given by Eq. (39). Therefore, as shown 
in Fig. 4(b), a typical iso- value surface of C(k) is the 
superposition of two crossed wheel-shaped spectra. 

2. Fully coherent case 

In the coherent-speckles case the 3D-correlation func- 
tion in real space and the corresponding power spectrum 
C(k) are the same as for the incoherent-speckles case, 
plus a coherence term. We then have 

C(k) = (V R /2) 2 [ci sp (k x ,k y , k z ) + c lsp (k z ,k y , k x ) 

~r~2c C oh (Me i ^yi kg)} (41) 

where 5i sp (k) is given by Eq. (39), c co h(k) is the Fourier 
transform of 



Ccoh(r) = Jci sp (x,y,z) x a sp (z,y,x) 

(l + 4ff)cos[0(r)]+2^sh#(r)] 

^1 + 4^/^^1 + 4x2/^ 

and <t>(r) = %f-(x - z) - T^rj^ ~ 7^ i+4z 2 %j| ' 
where ci sp (r) is given by Eq. (36) with w x — w y = w. 
The latter term mainly creates two broad structures 
(bumps), centered on the kx = (k^ — k z )/\/2 axis [see 
Fig. 4(c)]. For the parameters of Ref. [58], (T»/a± ~ 5.8 
and A L /cr^ ~ 2.16, these bumps are located at kx — 
+3.8CT7 1 . 



IV. SINGLE-SCATTERING 

We now focus on the first time scale introduced in 
Sec. II A: The scattering mean free time. 

A. Scattering mean-free time 

In order to calculate the scattering mean free time, 
defined previously in Eq. (14), we retain only the low- 
est order contribution to the self-energy (Born approxi- 
mation). Within this approximation, the Born series of 
Eqs. (5)- (6) is truncated after the first two terms, which, 
according to Eq. (8), yields 

£(£) = VG (E)V. (42) 

For homogeneous disorder, (k|E(JS)|k') = (27r) d (5(k — 
k')£(£,k) with 

f dk" - 

E(£7,k) = J C(k-k")G (£,k"), (43) 

where C(k) is the disorder power spectrum (Fourier 
transform of the correlation function [95]). Using 
Eq. (14) and the disorder- free Green function, Eq. (9), 
we thus have 

hi , 

rJS.k) = — -= — , 44 

2^<C(k-k')) k , |£ ' 

where 

("-W = /^ a -*[^-^)] («) 

represents the integration over the k-space shell defined 
by e(k) = E. Equation (44) allows one to determine the 
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Figure 5. (Color online) On-shell scattering mean free time 
Tg £ = r s (E, k E k) [Eq. (46) for |kj = k E ] along the k x 

(solid red line) and k H directions (dotted blue line) for the 
2D speckle potential defined in Sec. Ill A with £ = 4. The 
solid black lines are the isotropic low-energy limit obtained 
for k E a± <C 1 [Eq. (47)] and the high-energy limit obtained 
for k E a± 3> £ [Eq. (48)]. The insets show the angular depen- 
dance of t e j at two different energies [with the parametriza- 

tion k = (cos 8, sin 6)]. The points on the lines are color- and 
shape-coded to match those in the insets. 



scattering time from the two-point correlation function 
of the disorder. In the following we discuss anisotropic 
properties of the scattering time for the 2D case (the 3D 
cases are presented in Sec. VII A). 

In the case of isotropic disorder [i.e. C(k— k') = C(|k— 
k'|)] the scattering time does not depend on the direction 
of the incoming wave vector k. In general, the scattering 
is however anisotropic: (7(|k — k'|) = C(\k — k"|) and the 
probability that the particle acquires a direction k' or k" 
after the single-scattering event are different. Isotropic 
scattering is found for <5-correlated disorder, which is a 
key model of disordered metals. In this case, each scatter- 
ing event leads to the complete loss of the initial momen- 
tum direction k. In the case of anisotropic disorder we 
are interested in, not only the scattering is anisotropic, 
but it also depends on the direction of the incoming wave 
k. 



where k = k/|k| is the unit vector pointing in the direc- 
tion of k, ft^ is the k-space solid angle, k E = \j2mE/h is 
the momentum associated to energy E in free space and 
E a± = h 2 jma\ is the correlation energy of the disorder. 
The scattering time (46) is plotted in Fig. 5 as a function 
of energy along the two main axes, for |k| = k E and for a 
fixed geometrical anisotropy £ = 4. We use the notation 



•em 



T s (E,k E k.). Let us discuss some limiting cases. 



In the low-energy limit, k E a± -C 1, we have 



KE„ 



'EM 







TT 

E 2 



HE 




eE 2 



(47) 



which is displayed in Fig. 5 (left-hand side black lines). 
In this limit the de Broglie wavelength of the parti- 
cle (2w/k E ) exceeds the correlation lengths of the dis- 
order (cr± x and ua_ v ) and the speckle can be approxi- 
mated by a white-noise (uncorrelated) disordered poten- 

~ 2 

tial. More precisely, Eq. (37) becomes C(k) ~ V 2 tt^t 
(see Sec. Ill A) and t e f. is isotropic, constant, and it only 
depends on the product V 2 <T± x cr± y (up to corrections of 
relative order E/E a± ). 

In the opposite, high-energy limit, k E a± 3> £, the de 
Broglie wavelength of the particle is much smaller than 
the smallest correlation length of the disorder. The par- 
ticle then behaves 'classically'. Since C*(k) has a wider 
extension in the k y direction than in the k^ direction (for 
£ > 1), there are more scattering channels for particles 



travelling along x so that t e ^ 
we find 



< T 



EMv 



More precisely, 



'EM 



frE a± k E a ± 
V 2 r ^ 



U2 _i_ C2U2 
K x + S K v ' 



(48) 



which is shown in Fig. 5 (right-hand side black lines) . In 
particular, we find that in the high-energy limit r 



EM 



oc 



E. 



It is also interesting to study the anisotropy factor of 
the scattering time 



6 



'BM, 



EMv 



(49) 



B. Anisotropic Gaussian speckle (2D) 

Let us consider the 2D anisotropic speckle poten- 
tial of geometrical anisotropy factor £ introduced in 
Sec. Ill A. Replacing C(k) by Eq. (37) in Eq. (44) and 
using the disorder-free dispersion relation of the vacuum 
in Eq. (45), we obtain the scattering mean free time 



hE„ 



2t 



which is shown in Fig. 6 as a function of E/E a± and £. 
As already mentioned t e ^ is isotropic in the white-noise 
limit, so that £ s ~ 1 for k E a± <C 1 (left-hand side red line 
in Fig. 6). When increasing the energy, the scattering 
time first increases along the direction with the largest 
correlation length, i.e. the direction in which C(k) is nar- 
rower (x for £ > 1). Therefore, £ s increases with E, for 
sufficiently small values of E/E a± , and we have £ s > 1. 
Using Eq. (47), an explicit calculation yields 
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Figure 6. (Color online) Anisotropic factor £ s = t e j /t e g 
as a function of E/E a± and 5, for the 2D speckle potential of 
Sec. Ill A. The red lines are the low (£ s — s> 1) and high energy 
limits (£ s -> |) [see Eqs. (47) and (51)]. 



For k E <Tj_ ^> £, using Eq. (48), we obtain 



& - 7' ( 51 ) 




1 35 14 1 45 4 47 448 

k u a ± k u <j ± 



Figure 7. (Color online) On-shell spectral function as a 
function of k for the 2D speckle potential of Sec. Ill A, 
Vr = 0.2Ecr ± , and £ = 4. The top row shows the full spec- 
tral function. The bottom row shows cuts along the k x (solid 
red lines) and k y axis (dotted blue lines). The two columns 
refer to different energies: E — E a± (left) and E = 10E a± 
(right), which correspond to the dots and the squares in Fig 5, 
respectively. 



V. BOLTZMANN DIFFUSION 



which shows that the anisotropy factor of scattering is 
proportional to the inverse of the geometrical anisotropy 
(right-hand side red line in Fig. 6). Note that the clas- 
sical limit relation (51) is universal provided that the 
configuration anisotropy factor is well defined i.e. that 
the disorder correlation function can be obtained by the 
anisotropic homothety of an isotropic one, C(x,y) = 
Ciso(x,£,y)- In this high-energy limit, £ s < 1 (contrary 
to the low-energy limit case). Therefore, for any value of 
t e k exibits an inversion of anisotropy when the energy 
increases, typically at E ~ E a± . 

As described in section II C the scattering time is the 
width of the spectral function. It can be measured in a 
2D experiment such as that of Rcf. [53] by monitoring 
the momentum distribution of an almost energy-resolved 
wavepacket [66]. To illustrate this, a plot of the spectral 
function as a function of momentum and at fixed energy 
is shown in Fig. 7. In this plot the scattering time is 
calculated in the Born approximation (as above), and we 
have neglected the real part of the self-energy £'(£", k) in 
Eq. (13). In each direction k the spectral function peaks 
at At e k/fi and has a width proportional to 1/t e g.. The 
anisotropy of the scattering time is revealed in the angle- 
dependence of both these quantities. It is more apparent 
in the angular dependence of the amplitude, which shows 
marked peaks. At low energy, the maxima are located 
on the k x axis, while at high energy, they are located 
on the k y axis, which signals inversion of the scattering 
anisotropy. 



We now turn to the behaviour of the spatial density in 
the incoherent diffusive regime, which is characterized by 
the Boltzmann diffusion tensor Db(-E). We first give an 
explicit formula for the diffusion tensor, in the framework 
of the usual on-shell approximation, and then apply it to 
2D disorder (3D cases are discussed in Sec. VII B). 



A. Solution of the Bethe-Salpeter equation 

In the independent scattering (Boltzmann) and weak 
disorder (Born) approximation, only the first term in 
Eq. (22) is retained and the irreducible vertex function 
U equals the disorder structure factor [81]: U ~ Ub = 
V ®V and 

U k , k > (q, oj, E) ~ U BKW = <5(k - k'), (52) 
or equivalently 



U B 



i 
i 

1 



(53) 



Then, incorporating Eq. (52)-(53) into the BSE (20)-(21) 
and expanding it in series of U, one finds 



k' 



<i> 



k 



k' 



— P^ 


r 


-Pi— 






-^— 



(54) 
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where the diffuson T reduces to ladder diagrams: 



• •■ 

I I 

1+ 1 

I I 

i 



' + 



(55) 



It describes an infinite series of independent scattering 
events, which leads to Drudc-likc diffusion. 

In appendix A, explicit calculations are detailed. In 
brief, in the long-time (u — > 0) and large-distance (|q| — > 
0) limit the vertex $ is the sum of a regular term and a 
singular term [75, 98]: 

*k,k' (q, «, E) = $f k 8 , (q, u, E) + ^ (0, 0, E). (56) 



The regular part is given by 



*2t (o.o 



E)=Y — L- 
' 1-A' 



■fEM$EM. ( t>EM.'fEM', (57) 



where /^,k = /k(q = 0,w = 0, E 1 ) [see Eq. (25)] and 
k (^e) are tne eigenvectors (eigenvalues) of an integral 
operator involving the disorder correlation function and 

Sem [109]: 



dk' 

(2tt) 



5 C(k-k')/ £ , k ' #|,k' 



a e9em- 



(58) 



The regular part contributes to the finite time and finite 
distance propagation of the density which we disregard 
here. The singular part is more interesting. Its existence 
is a direct consequence of the Ward identity [110] which 
expresses the conservation of particle number, and which 
guarantees that one of the eigenvalues of Eq. (58) is equal 
to one A^ =1 = 1 (see appendix A). In the framework of 
the on-shell approximation, such that e(k) = e(k ) = E, 
in the long time and large distance limit (|q|, u>) — > 0, the 
vertex $ is given by 



2ir 7k(q,£) 7k '(q,£) 



HN (E) -?:cj + q-D B (£)-q 
with Nq(E) the disorder-free density of states, and 



7k(q,£) = 



A (E,k) 



2tt 



1 - 



2ni 



(59) 



(60) 



E,i.'^E.k'' k '\ E 



where Ao(E, k) = 2ir5[E — e(k)] is the disorder-free spec- 
tral function. Equation (59) shows that the vertex $ 
is dominated by the diffusion pole (ihui — hq-T>B(E) ■ 
q) _1 . The Boltzmann diffusion tensor Db(£) has com- 
ponents [75] 



D^(E) 



N (E) 

Y ^ T-a 



k|_E 



(61) 



2?r ^ A| 



EM 



k|E 



£,k^^, k 



hki/r 



•em 



T s (E,k E k) = h/2Tr(C(k E k 



where i 

k'))k'|B is the on-shell scattering mean free time [see 
Eq. (44)], and (...)k|_B represents integration over the k- 
space shell defined by e(k) = E [see Eq. (45)]. The func- 
tions 0™ r and the real- valued positive numbers A^ are 

the solutions of the integral eigenproblem (58), which be- 
comes, in the on-shell approximation (sec appendix A), 



^(r^C(a-k')^) k , 



n jji 
EM 



(62) 



normalized by f < T £,k^,k^ k ) k |B = S n,m [75]. It fol- 
lows from Eq. (61) that the incoherent (Boltzmann) dif- 
fusion tensor Db (E) is obtained from the two-point dis- 
order correlation function C (r) , which determines t e £ 
[see Eq. (44)] as well as 0™ - and A| [see Eq. (62)]. 

In the isotropic case (for details see appendix B), 
Eq. (62) is solved by the cylindrical (2D) or spherical 
(3D) harmonics Y, m , the same level harmonics [i.e. with 
the same I] being degenerate in A^. Then, it follows 
from the symmetries of the cylindrical/spherical harmon- 
ics that only the first term in Eq. (61) plus the p-level 
harmonics {Zf 1 in 2D and Yf 1 with m = -1, 0, 1 in 3D; 
the cylindrical harmonics Zf~ x are defined in appendix B) 
couple to v and contribute to D B (E). Incorporating the 
explicit formulas for <^>™ ^ and X E [see Eqs. (Bl) to (B7)], 
we then recover well-known expressions for isotropic dis- 
order [61-63, 65]. 

In the anisotropic case, harmonics couple, and the <\> n ^ 

are no longer cylindrical/spherical harmonics. Then the 
calculation of the diffusion tensor requires first the diag- 
onalization of the integral operator (62) whose solutions 
are then incorporated in Eq. (61). In the following, this 
is done numerically for an anisotropy factor relevant to 
current matter- wave experiments. 



B. Anisotropic Gaussian speckle (2D) 

Consider again the 2D anisotropic speckle potential of 
Sec. Ill A. The first step in the calculation of Db is to de- 
termine the eigenfunctions 6"l - and the associated eigen- 

T em 

values A^ of Eq. (62). We solve Eq. (62) numerically, by 
a standard algorithm of diagonalization, with 2 9 = 512 
points, regularly spaced on the k-spacc shell |k| = k E . 
The diffusion tensor is diagonal in the basis made by the 
symmetry axes of the correlation function (37): {u x , u y }. 

The eigenvalues and some eigenfunctions obtained nu- 
merically are shown in Fig. 8 for various values of E/E a± . 
As discussed above, we find A^ =1 = 1. For E <C E a± , 
only the first term in the right-hand side of Eq. (61) con- 



tributes to the diffusion tensor since all A^ 1 are vanish- 
ingly small. When the energy increases, the values of the 
coefficients A^ 1 increase. It corresponds to an increase 
of the weight of the terms associated to the orbitals with 
n > 1 in Eq. (61), and a priori all the orbitals with 
n > 1 might have an increasing contribution. However, 
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Figure 8. (Color online) Top row: Eigenvalues of Eq. (62). 
Bottom row: Angular dependence of the eigenfunctions 0™ ^ 
for n = 1 (dashed black line), 2 (solid red line) and 3 (dotted 
blue line). We use the parametrization k = (cos 9, sin 0) for 
the 2D speckle potential of Sec. Ill A with £ = 4. The different 
columns refer to different energies (indicated on top of the 
figure) . 



we find that, the symmetry properties of the functions 
d> n - cancel the contributions of most of them, and only 

the orbitals with n = 2 and 3 do contribute (see below). 

In the low energy limit, one can develop Eq. (37) in 
powers of |k|. Up to order 0(E 2 /£ 4 E 2 ± ), the first three 
eigenfunctions are given by: 




Figure 9. (Color online) Components of the diffusion ten- 
sor: (soild red line) and (dotted blue line) for the 
2D speckle potential of Sec. Ill A and £ = 4. Solid black 
lines are limit values at small E/Ecr ± [Eqs. (66) and (67)], 
with the isotropic white-noise limit Dg(E) = D^(E) ~ 
fr£,EE a± /m%V£. For large E/E a± we find D B (E) ~ E 5/2 
(see text); a fit of the numerical data gives the prefactors 
Dl = 4.43£ 5/2 /% 2 £^ 2 and Dl = 1.24 E b/2 /VIeI 1 * (see 
dotted black lines). The inset shows the transport anisotropy 
factor £ B = Dl/Dl. 
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2i?E a 
with eigenvalue A^ = 1; 



i + (£ 2 - 



V2 + B 2 
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E- 



C 4 E 2 



(63) 
(64) 



with eigenvalue \ 2 E = E/2E a± ; and 



V2 + Bi 



E 







E 2 



(65) 



with eigenvalue X% = E/2£ 2 E a± , where B2 and B3 
are constant values that do not intervene in the follow- 
ing. In this limit the numerical results agree very well 
with the analytical findings (which for clarity are not 
shown on Fig. 8). In the very low energy limit, the dis- 
order power spectrum becomes isotropic and constant, 
C*(k) ~ V 2 tt cri/£, [see Sec. Ill A and Eq. (37)]. The 
orbitals <f>" ~ are thus proportional to the cylindrical 

harmonics, which are exact solutions of Eq. (62) in the 
isotropic case (see appendix B, and use the parametriza- 
tion k x = cos 9 and k y = sin#). In contrast to the 
isotropic case where the values of \ E are degenerated 
in a given /-level, here we find that the degeneracy inside 
a I level is lifted for any anisotropy £ 7^ 1 [see the val- 
ues of A^ below Eqs. (64) and (65)]. When the energy 
further increases, the anisotropy plays a more important 
role and the harmonics are more and more distorted (see 
Fig. 8). However their topology remains the same, and 



in particular the number of nodal points and their posi- 
tions are unchanged. In the following, we thus refer to 



^, ±1 -like orbitals. 



E ' 4>eI and T EM U1 



Incorporating the values of A 

Eq. (61), we can determine the Boltzmann diffusion ten- 
sor. Figure 9 shows the resulting eigencomponents of the 
diffusion tensor. In the low energy limit (E <C E a± ), us- 
ing Eqs. (47), (64) and (65), we find that the first term 
in the right-hand side of Eq. (61) gives the leading con- 
tribution to Db(£) (of order E/E^^. This contribution 
is isotropic owing to the isotropy of r E ? at low energy 
and of the underlying medium. At very low energy, in 
the white-noise limit, we recover an isotropic diffusion 
tensor D%{E) = Dl(E) ~ h£EE a± /rrnrV 2 . The scaling 
D'£(E) oc E is universal for 2D disorder in the white- 
noise limit (when it exists). The -like orbital 4> E -^ 

contributes to the next order of D* and the Z^ 1 -\ike or- 
bital cj) 3 E fc to Dl. Up to order 0(E 3 /£ 6 E* ± ), we obtain 



DUE) 



and 



DUE) = 



771 V 2 



hE^ 

m V 2 



jE 

nEcr. 



eg 

nE ni 



E 2 9£ 2 + 3 



8£ 



E 2 3£ 2 



*E 2 ± 



8£ 











E : '' 



i & El^ 
(66) 



E :i 



A & El 
(67) 

which are displayed on Fig. 9 (left-hand side solid lines). 
When the energy increases, the anisotropy first comes 
from the anisotropic contribution of the scattering time 
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Figure 10. (Color online) Boltzmann transport anisotropy 
factor £b = D^/D^ as a function of E/E a± and £ for the 
2D speckle potential of Sec. Ill A. The inset shows the high 
energy asymptotic value (cut at E = lQ*E a± ). The dotted 
red line in both the figure and the inset is £. 



t e k' an< ^ from the lift of the degeneracy between 
and A|;. When the energy further increases, the harmon- 
ics are distorted, - but their symmetries (i.e. periodicity 
and parity) are preserved (see Fig. 8). Hence, for the 
same reasons as in the isotropic case (see appendix B) 
only the Z^-like orbitals couple to v in Eq. (61) and 
contribute to D B while the others don't. The associated 
increase (see Fig. 8), the weight of the second term in 
Eq. (61) increases, and the components of the diffusion 
tensor show a very different behavior in the large- E limit. 
For k E a± S> £, we found t e £ oc k E (see Sec. IV B). In 
addition, we find numerically a weak topological change 
of the orbitals with energy for E/E CT± > 10 2 . Therefore 
the evaluation of Db with E is mainly determined by 
the normalization condition [see formula below Eq. (61)], 
which yields 0™ ^ oc 1/ yku. Then, assuming the scal- 
ing 1 — oc 1/E, also verified numerically, we obtain 
D^{E) oc _E 5 / 2 , which matches the numerical results (see 
dotted black lines in Fig. 9). This scaling is similar 
to that found for isotropic disorder [62]. As shown in 
Fig. 9, the change of slope between the low- and high- 
energy regimes is different in the two directions. For 
this reason, the anisotropy factor of the diffusion tensor, 
£b = D*/Db shows a nonmonotonous behaviour versus 
E, with a marked peak (see inset of Fig. 9). 

The Boltzmann transport anisotropy factor £b is 
shown in Fig. 10 for various geometrical anisotropics £. 
As it is well-known, the scattering and transport mean 
free times are different quantities in correlated disorder, 
due to angle-dependent scattering [81, 111, 112]. In par- 
ticular, in the 2D speckle we consider, we do not find 
any inversion of the anisotropy of the diffusion, con- 
trary to the scattering time, i.e. the component D*(E) 
of the diffusion tensor is always larger than the compo- 
nent Db(E). For large values of E/E a± , the Boltzmann 



transport anisotropy £b reaches a constant value (see the 
inset of Fig. 9 for a cut at £ = 4), which increases with 
the geometrical anisotropy £ (see inset of Fig. 10). This 
asymptotic value is larger than £ for small £ and smaller 
for larger £. Therefore the anisotropy of the diffusion in 
the classical regime is not simply related to the spatial 
anisotropy. 

The two distinct regimes found in the behaviour of Db 
and the non-trivial anisotropy effects make the Boltz- 
mann diffusion regime in anisotropic 2D potentials very 
interesting for future experiments. Those properties 
could be probed by imaging directly the atoms in the 
2D speckle (as in ref. [53]) and controlling the width of 
the atomic energy distribution. 

VI. WEAK AND STRONG LOCALIZATION 

Having discussed the incoherent (Boltzmann) trans- 
port properties, we now consider interference effects, 
which lead to weak and strong localization. We first 
describe the quantum corrections (Sec. VIA), then the 
self-consistent theory (Sec. VI B), and apply it to the 2D 
speckle potential (Sec. VIC). The 3D cases, which follow 
the same route, are discussed in Sec. VII C. 

A. Weak localization correction 

We calculate corrections to Boltzmann diffusion by 
taking into account quantum interference terms between 
the multiple-scattering paths. Those interferences ap- 
pear when the correlated scattering events do not occur 
in the same order in the propagation of the field and 
its conjuguate. This is diagrammatically translated into 
crossing correlation lines as in the second term of Eq. (22) 
for example. In the weak scattering regime only the two- 
point correlations are retained in the scattering diagrams 
and the leading scale-dependent corrections to the clas- 
sical conductivity are given by the maximally crossed di- 
agrams [6, 75, 98, 113]: the cooperon [Eq. (68)] and the 
first two Hikami boxes [Eqs. (69) and (70)]. 
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where the cooperon X is the sum of maximally crossed 
diagrams 




\ / 
\ / 

A 
/ \ 



I 



\ / 



and 




(71) 



(72) 



is the renormalized vertex function (see appendix C 2). 

Using time-reversal invariance [23, 81, 110, 114, 115], 
the cooperon X can be expressed in terms of the diffuson 
r [defined in Eq. (55)] 



X ky {q,u J ,E)=T l 



_ a (k + k',w,£). (73) 



The diffusion pole carried by T in the limit (w,q) — > 
leads to a divergence of X when to, k + k' — > 0. In 
appendix C 3 we translate those diagrams into equations, 
and show that 



Act(uj,E) 



dQ 



1 



wN {E) J (2ir) d -ihuj + HQ ■ D B (£) • Q' 

(74) 

Using Einstein's relation (34) we then obtain the dy- 
namic diffusion tensor D»(w,£) = Db(£) + AT>(uj,E), 
with [75] 



£D(u),E) 
D B (£) 



dQ 



1 



ttNq(E) J (2ir) d -ihoj + HQ ■ T) B (E) ■ Q' 

(75) 

Note that the quantum corrections AD(cj, E) do not ex- 
plicitly depend on the disorder [i.e. on C(k)], but only 
on the Boltzmann diffusion tensor Db(-E) [75]. In other 
words, in this approach, Boltzmann incoherent diffusion 
sets a diffusing medium, which contains all necessary in- 
formation to compute coherent terms [116]. In particu- 
lar, it follows from Eq. (75) that the weak localization 
quantum correction tensor AD(w, E) has the same eige- 
naxes and anisotropics as the Boltzmann diffusion ten- 
sor Db(£). Thus the anisotropy can be removed by 
rcscaling distances along the transport eigenaxes u by 
J I D^ v (i.e. momenta are rescaled by y/D™ /D™) with 
D| v = det{DB} 1 ^ d the geometric average of the Boltz- 
mann diffusion constants. Since AD is always negative 
in the limit u> — > + , the weak localization correction 
determines slower diffusion than the one obtained from 
incoherent diffusion. Equivalently, as long as the correc- 
tion (75) is small, one can write 



D B (£) 



= 1- 



1 



dQ 



1 



n*(u,E) ttNq(E) J (2n) d -ihaj + HQ -T) B (E) ■ Q' 

(76) 

which is the lowest-order term of a perturbative expan- 
sion of 1/D»(w, E). 



B. Strong localization 

The quantum interference correction (75) has been de- 
rived perturbativcly and is therefore valid as long as the 
correction itself is small, i.e. for T)-q(E) — D*(lj,E) <C 
Db(£). In order to extend this approach and eventually 
describe the localization regime where D* vanishes, Voll- 
hardt and Wolfle [110, 114] proposed to self-consistently 
replace Db(E) by the dynamic diffusion tensor D*(w, E) 
in the right-hand side of Eq. (76). For isotropic scatter- 
ing this procedure amounts to resumming more divergent 
diagrams than the cooperon (which contain a square of a 
diffusion pole), thus contributing to localization [80, 110]. 
Generalizing this standard approach to anisotropic dis- 
order yields 



Dt 



= 1 



1 



dQ 



1 



ttN (E) J (2ir) d -ihuj + HQ ■ D*(u;) • Q' 

. (77) 

In dimension d > 2 the integral in the right-hand side of 
Eq. (77) features ultraviolet divergence, i.e. in the limit 
|Q| — > +oo. Since the diffusive dynamics is relevant 
only on length scales larger than the Boltzmann mean 
free path Vq{E) = dy / m/2E D^(E) along each transport 
eigenaxis, we regularize this divergence by setting an up- 
per ellipsoidal cut-off of radii 1 /Zg in the integral domain. 
It corresponds to an isotropic cut-off in the space rescaled 
according to the anisotropy factors of D B as described 
above. 



C. Anisotropic Gaussian speckle (2D) 

We now solve the self-consistent equation (77) for the 
inverse dynamic diffusion tensor in the 2D case [110, 114]. 
In the long time limit u> — > + , the unique solution of 
Eq. (77) is of the form D,(w,S) - 0+ -iuLf oc {E), where 
Lioc(-S) is a real positive definite tensor. As described in 
Sec. II D, it leads to the exponentially decreasing prop- 
agation kernel (30). Solving Eq. (77) then yields the 
anisotropic localization tensor, 




1/2 



(78) 



where Z| v (£) = d^fm/2ED™{E). The eigenaxes of the 
localization tensor are thus the same as those of the 
Boltzmann diffusion tensor and its anisotropy factors are 
the square root of those of D B (E), i.e. £i oc = Lf /£f oc = 

We now apply the self-consistent theory to our run- 
ning example: the 2D anisotropic speckle potential with 
correlation function given by Eq. (37). Including the re- 
sults for the Boltzmann diffusion tensor Db(-E) obtained 
in Sec. VB into Eq. (78), we find the localization tensor 
Lioc (■£')■ Figure 11 presents the eigencomponents of Li oc 
in its eigenbasis {u x ,iiy} as a function of energy, for a 
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Eq. (80), the 2D localization length increases exponen- 
tially at large energy (hence the limited energy range in 
Fig. 11). Therefore measuring it experimentally with ul- 
tracold atoms [54, 117, 118] is very challenging and can 
be done in a rather narrow energy window, in which i loc 
is larger than the resolution of the imaging system (L Tes ) 
but smaller than the size of the sample (L sys ). This is il- 
lustrated for ct_l = 0.25 fim on Fig. 11 by the grey dashed 
lines L rcs ~ 15 /xm and L sys ~ 2 mm, which are typical 
values extracted from Refs. [53, 58]. 



VII. THREE-DIMENSIONAL DISORDER WITH 
STRUCTURED CORRELATIONS 



Figure 11. (Color online) Components of the localization ten- 
sor Lf oc (solid red line) and Lf oc (dotted blue line) for the 2D 
speckle potential of Sec. Ill A, with £ = 4 and V R = 0.2E a± 
and 2E a± . The solid black lines are the limiting behaviour for 
small values of EjE^^ [Eq. (79)] and the dotted ones for high 
values of E/E a± [Eq. (80)]. The dashed grey lines indicate 
typical values of the imaging resolution (L res ) and the system 
size (L aya ) in ultracold-atom experiments, see text at the end 
of Sec. VIC. 



geometrical anisotropy of £ = 4 and two different ampli- 
tudes of the disorder V R /E <7± =0.2 and 2. At low energy 
(E < E (J± ,V R ,V 2 /E <7± ), using Eqs. (66) and (67), we 
find 



E 3 f 3 / 2 2E 
L^{E)=a^^ 



V 3 it E a 



1 



jEE a 
2V 2 



E (18±3)£ 2 + (18t3) 



E, 
O 



E 2 



16£ 2 
E 2 



EE 2 



(79) 



where the upper sign holds for direction x, and the 
lower sign for direction y. Equation (79) corresponds 
to the solid black lines in Fig. 11. As Db is almost 
isotropic for E/E a± < 1 (see Fig. 9), Li oc is also al- 
most isotropic in the whole range presented in Fig. 11. 
Equation (79) describes an isotropic localization tensor 
with an anisotropic correction which is significant only 

if V R /E a± > £ 3/2 /VF--l (- 2 for f = «• At hi g her 
energy, when k E l^(E) = 2mD™(E)/h > 1, we expect 



LISE) ~ a x 



2 my/D™(E)D2(E) 
k K a\ h 



,mnD B "(E)/h 



(80) 

which is plotted as dotted black lines in Fig. 9. Ac- 
cording to Eqs. (66) and (67) (retaining only the 
lowest-energy term), this regime appears for E/E a± > 
(Tr/2£)(V R /E ry± ) 2 . When £ = 4 (as in Fig. 11), it gives 
E/E a± > 0.015 for V R /E a± = 0.2 and E/E a± > 1.5 
for V R /E cr± = 2. As predicted by the scaling theory 
of Anderson Localization [23] and explicitely seen in 



In this section we apply the formalism introduced 
in Sees. IV to VI to the 3D speckle potentials 
of Sees. IIIB and IIIC. We discuss single-scattering 
(Sec. VII A), Boltzmann diffusion (Sec. VII B) and local- 
ization (Sec. VII C) properties, successively for the single- 
speckle and orthogonally-crossed-speckles configurations. 
We recall that those systems are relevant for ultracold 
atoms experiments. In particular our configurations ap- 
ply to Ref. [57] (single-speckle) and Rcf. [58] (coherent 
orthogonally-crossed speckles), respectively. This section 
can be viewed as a detailed version of Ref. [79] . 



A. Single-scattering 

1. Single speckle configuration (3D) 

Let us first consider the single- speckle case. Inserting 
Eqs. (38) and (39) into Eq. (44), we find the scattering 
mean free time 



r s (£,k) 



(27r) 2 /k E a ± 



V, 



R J df2 fc ' c lsp (k E k - k)/o^ 



(81) 



which is shown in Fig. 12 for |k| = k E [as for the 2D case 
of Sec. IV B we define the on-shell scattering mean free 
time t e f = t s (E, k E \i)]. Since C(k) is isotropic in the 
(k x ,k y ) plane, t e g. only depends on the polar angle 6 

between k and k z and not on the azimutal angle 4>. We 
find that the scattering time is an increasing function of 
energy. It is also shorter for particles travelling along the 
z direction (t e g < t e £ ) for all values of E. As 
for the 2D case analyzed in Sec. IV B, this is due to the 
wider extension of C(k) in the plane (k x , k y ), which offers 
more scattering channels to particles travelling along z. 
In contrast to the 2D speckle case however, t e g shows 
no inversion of anisotropy. 

In the low energy limit (k E a± <C 1), t e g converges to 
a constant value. In contrast to the 2D case, it signals 
the absence of a 3D white- noise limit [119]. This can be 
attributed to the strong anisotropic divergence of C(k) 
when |k| — > which reflects the long-range correlations of 
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Figure 12. (Color online) Scattering mean free time t e £ in 
the 3D single- speckle case [Eq. (81)] with |k| = k E for the 
parameters of Fig. 4, in the (kz,k v ) plane (solid red line) and 
along the k z direction (dotted blue line). The black lines are 
the low-energy [k E a± <C 1, see Eq. (83)] and the high energy 
[k E a± 2> 1, see Eq. (84)] limits. Note that in both limits t e ^ 
is anisotropic, although for k E a± -C 1, the anisotropy is very 
small, £ s ~ 1.002. The insets show the angular dependence 
of t e £ at different energies [with the parametrization k = 

(ka, , k y , k z ) = (sin 6 cos <j>, sin 6 sin <j>, cos 6)]. The points on 
the lines are color- and shape-coded to match those in the 
insets. 



the disorder (see Sec. MB). More precisely, for |k|erj_ <C 
1, we have 



gl ( k ) ~ it 3 / 2 ^a-(k) - n 3/2^i^ 



Replacing c isp in Eq. (81) we then find 



hE a± 4v^F 



4 W/ W+R, 



£•2 _i_ £.2 

K x ' K y 



r J dn u - c(k - k) ' 



(82) 



(83) 



which is independent of E. Equation (83) is plotted as 
solid black lines on the left-hand side of Fig. 12. Note 
that t e g does not become strictly isotropic in this limit. 
However, the residual anisotropy of the scattering time, 
found from Eq. (83) and from the anisotropy of c(k) 
in Eq. (82), is very small, and practically unobservable 
( t e k { } / r £k — 1-002). When the energy increases, 
the scattering time in the (x, y) plane is the first to 
deviate significantly from the low-energy behaviour at 
E ~ E a (= 3 x 10~ 2 E a± for the parameters of Fig. 12), 
while the scattering time in the z direction increases only 
at E ~ E a± . This can be understood again by the nar- 
rower width of the power spectrum C (k) in the k z direc- 
tion. 

In the high-energy limit (k E a± 3> 1) the k-space shell 
integral of Eq. (81), which is done on a sphere of radius 



k E containing the origin, can be reduced to integrating 
c iS p on the plane which is tangent to the sphere at the 
origin. We then find 



A^Kk E (jA 



'em 



(84) 

In particular, we find t e = hE a± k E cr±/2V 2 y/ir, 

hE ai k E o\lV 2 ito\\ (both shown as the right- 



hand side solid black lines in Fig. 12). The anisotropy of 
the scattering then becomes significant for the parame- 



ters of Fig. 12, r 



E.k 



e i ~ v^" "!! / 2ctj^ in this limit. 
The high-energy scaling t e &oc k E , which was also found 
in our 2D speckle, is quite universal: as long as the power 
spectrum is of finite integral in all the planes (lines in 2D) 
crossing the origin, the procedure described above can be 
applied to Eq. (44). Then t e £ only depends on the dis- 
persion relation e(k) and, in particular, it is independent 
of the space dimension. 



2. Orthogonally-crossed speckles (3D) 

We now consider the case of two orthogonally crossed 
speckle fields, that can be either mutually incoherent 
or coherent, and whose power spectrum are given by 
Eqs. (40) and (41) respectively. The 3D on-shell scatter- 
ing mean free time [Eq. (44) with |k| = k E ] is presented 
in Fig. 13 in both configurations. 

The power spectrum of the incoherent- speckles case 
is made of two orthogonally-oriented spectra, similar to 
that of the single-speckle case. As a consequence, its scat- 
tering time [see Fig. 13(a)] is qualitatively similar to the 
single-speckle one. It shows two disctinct regimes: t e £ 

constant at low energy and t e g. oc \f~E~ at high energy. 
Note that even though the directions X and Z are equiv- 
alent, t e g has a dependence in 6 = (k, kz) and in ip (the 
azimutal angle in the (kx, ky) plane) because the corre- 
lation fonction does not show rotation invariance around 
any axis. As an example t e g is also presented on 

Fig. 13(a). Note also that the anisotropy between the di- 
rections of minimal (kr x z y) and maximal (ky) scattering 
times is reduced compared to the single-speckle case. 

In the coherent- speckles case, the power spectrum 
is the incoherent one plus a coherence term that cre- 
ates two additional bumps centered on the kx axis [at 
kx — ±3.8er^ 1 for our parameters, see Fig. 4], As al- 
ready mentionned, the scattering time (taken on-shell) 
in the Born approximation samples C(k) on a k-space 
shell of radius k E centered in fc B k [see Eq. (44)]. There- 
fore we recover the incoherent-speckles case at low en- 
ergy. The bumps play a role for 2k E a± > 3.8, i.e. for 
E > 1.8E a± . They offer more scattering channels to 
the particle travelling along Y and Z, making the cor- 
responding scattering times drop below the values ob- 
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Figure 14. (Color online) Single-speckle case. Eigenvalues of 
Eq. (62) at various energies indicated on the figure (top row). 
Topography of the eigenvectors </>™ £, at the same energies, 

which mainly contribute to D^. (bottom row), D B (2 nd row) 
and Dg (3 rd row) respectively [with the parametrization k = 
(k x ,\t y ,]i z ) = (sin 8 cos cj>, sin 6 sin <j>, cos 9)] . The values of n 
are indicated on the figure, the red lines locate the nodal 
lines. The points are color- and shape-coded to match those 
of Fig. 12. 



Figure 13. (Color online) Scattering mean free time t e £ in 
the 3D incoherent (top) and coherent- (bottom) speckles cases 
for the parameters of Fig. 4, along the symmetry axes of 
the correlation functions (see Fig. 4), kx = (k x — k z )/\/2 
(solid red line), ky = k^ (dotted blue line), and kz = 
(k. x + ]i z )/^/2 (dashed gray line) directions. The solid black 
lines are the low-energy limits obtained for k E a± <C 1 and 
the high energy limits obtained for k E a± ^> 1 [see Eq. (85)]. 
The insets show the angular dependance of t e £ at differ- 
ent energies [with the parametrization k = (kx,ky,kz) = 
(sin 6 cos <p, sin <d sin <p, cos 0)]. Note that this parametriza- 
tion differs from that of Fig. 12 because the symmetry axes 
of C(k) are different. The points on the lines are color- and 
shape-coded to match those in the insets. 



taincd in the coherent case [see Fig. 13(b)]. This leads 
to an inversion of the anisotropy of the scattering time 
with energy: t e ^ x < t E £ y ,t b £ z at low energy and 
t em x > t em y » T E,k z at hi § h energy. 

At low energy (k E a± <C 1) the absence of white-noise 
limit and the scaling of Ci sp (k) presented in Eq. (82) 
gives, as for the single-speckle configuration, a constant 
slightly anisotropic scattering time in both cases (the 
term c co h(k) present in the coherent case being negli- 
gible), with a very small anisotropy. When E increases, 
the scattering time in all directions deviates from the 
low-energy behaviour around E ~ E a (= 3 x lQ~ 2 E a± 
for the parameters of Fig. 13). In the high-energy regime 
(k E a± >• 1 for the incoherent and k E a± >• 3.8 for the 



coherent case) we have 

_ KE a± (2Tr) 2 k E <7 ± 

T E,k - y2 _ / ^ — \ 

R J cIkcIk' C I Kk z ,K', n^Jkl + fc2 J /V£a ± 

(85) 

which is displayed as the solid black lines on the right- 
hand side of Fig. 13. 

B. Boltzmann diffusion 

Let us now analyze Boltzmann diffusion in our 3D 
speckle potentials. It is obtained, as in the 2D case an- 
alyzed previously by solving Eq. (62) numerically and 
incorporating the results in Eq. (61). For the diagonal- 
ization of the integral operator (62) we use 2 7 x 2 7 = 
128 x 128 points regularly spaced on the k-space shell 
|k| = k E [120]. 

1. Single-speckle (3D) 

The eigenvalues X E of Eq. (62) for different energies, 
as well as the topography of the eigenvectors of Eq. (62) 
that dominate (bottom row), D\ (2 nd row), and 
(3 rd row) are shown in Fig. 14. We find (similarly as for 
the 2D speckle potential) that the </>" ^ are topologically 
similar to the spherical harmonics at all energies, i.e. they 
show similar nodal surfaces, but the associated X E are 
not degenerated in a given Z-like level. More precisely, 
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Figure 15. (Color online) Boltzmann diffusion coefficients 
along the transport eigenaxes (eigencomponents of Db) for 
the single-speckle configuration and for the parameters of 
Fig. 4. The dotted lines are power-law fits (D^ oc E lu ) to 
the data in the low and high energy limits. The inset show 
the transport anisotropy factor D^/D^,' y . 



due to the cylindrical symmetry of the power spectrum 
[see Fig. 4(a)], the value of X E associated to the V ; +m - 
like and Y" z -m -like orbitals are the same for a given m, 
but the degeneracy between the different values of \m\ is 
lifted. 

Figure 15 shows the resulting eigencomponents of the 
diffusion tensor in the single- speckle case. It is isotropic 
in the (x, y) plane, because of the rotation-invariance of 
the correlation function C(k) around the axis k 2 . For 
the same symmetry reasons as in the isotropic case (see 
appendix B) and as in the 2D case, only the p-lcvcl-likc 
orbitals couple to v. For k E a± <C 1, wc find that D% ,v is 
dominated by the first term in Eq. (61) and D* by the 
Y^-like orbital (n = 2 at all energies). For k E a± S> 1, the 
situation changes: while D* is still dominated by the Y®- 
like orbital, is now dominated by the Y 1 +1 -like orbitals 
and D\ by the Y 1 _1 -like orbitals (respectively n = 6 and 
5 at E = 50E a± in Fig. 14) with a contribution of the 
y 3 ±1 -like orbitals increasing with E [121]. Let us discuss 
the main features of Db ■ 

Firstly, as already discussed for the 2D case, the trans- 
port and scattering mean free times can be very different 
quantities in correlated disorder, and, in particular the 
anisotropy of Db can be very different from that of t e ^. 
Here we find that the diffusion tensor is larger along axis 
z (D* > Dl' v ) for all values of E (see Fig. 15), and the 
anisotropy of Db is thus reversed with respect to that of 
t e k ( we recau that we found t e ^ < t e ^ } for any 
E, see Sec. VII A). This is due to the fact that the (Yf- 
like) orbitals contributing to £>* are associated to values 
of X E larger than those contributing to D^' y (in Fig. 14, 
the 4> E £ are numbered by decreasing eigenvalues) 

Secondly, C(k) shows a strong anisotropic, infrared di- 
vergence in the paraxial approximation (see Sees. Ill B 



and VII A). Following-up with the scaling of ci sp (k), 
Eq. (82), used to show that E is independent of en- 
ergy for k E a± <C 1, and inserting it into Eq. (62) and the 
associated normalization, we find that A^ does not de- 
pend on E, and 0™ . is of the form ip n (k) / \fk~ E . Then, all 

terms in Eq. (61) are topologically unchanged and scale 
as E at low energy. The anisotropy of Db thus persists 
down to arbitrary low values of E and oc E, as ob- 
served in the left-hand side of Fig. 15 for k E a± <C 1 (i.e. 
E <C E a± ) . This is another manifestation of the absence 
of white- noise limit [122]. 

Thirdly, we found t e <■ oc \^E, and assuming weak 
topological change of the orbitals and the scaling 1 — A^ oc 
1/E (confirmed numerically), we get 4> E ^ l/k E and 

D™(E) oc E b / 2 . This scaling is confirmed in Fig. 15 by 
fits to the data for E ^> E a± (right-hand side dotted 
lines). This scaling was also found in our 2D example 
and for isotropic 3D speckle disorder (see Ref. [62] and 
appendix B). Remarkably, in spite of the different con- 
tributing terms in Eq. (61) at low and high values of E, 
the transport anisotropy is nearly independent of E with 
Dl/D*' y ~ 10 [see inset of Fig. 15]. 

2. Orthogonally-crossed speckles (3D) 

Let us turn to the crossed-speckles configurations, 
whose diffusion coefficients arc plotted in Fig. 16. Note 
first that in both the incoherent- and coherent-speckles 
configurations we recover the same general properties 
as for the single-speckle case, in particular the reversed 
anisotropies of scattering (r E ^ x z < t e £ y ) and diffu- 
sion (Dg ' Z > D\), the anisotropic suppression of the 
white-noise limit, and the scaling of the diffusion coef- 
ficients at low [D%(E) oc E] and high [D%(E) oc E 5 / 2 } 
energy (see left- and right-hand dotted lines in Fig. 16). 
Here however, the transport eigenaxes are the bisectors 
{X,Z} = (x=pz)/"\/2 and the axis Y = y, which are sym- 
metry axes for both correlation functions (see Fig. 4). 

In the incoherent- speckles case [Fig. 16(a)], Db is 
isotropic in the (X, Z) plane, even though the correlation 
function and the scattering time are not. This is due to 
the topology of C(k) which bears four symmetry axes in 
this plane (k-ck^kx and kz) [123]. We find that the 
0™ r contributing to the diffusion tensor are distorted 

compared to the single-speckle case but the number of 
nodal lines and their positions are still reminiscent of the 
17" spherical harmonics. In the end, the behaviour of 
the diffusion tensor is very similar to the single-speckle 
one and ~ > E>\. The transport anisotropy is 
nearly constant, but significantly reduced with respect to 
the single-speckle case, Dg' z jE)\ ~ 1.8. 

In the coherent- speckles configuration [Fig. 16(b)], as 
long as 2k E a± <C 3.8, the behavior of L>b(E) is gov- 
erned by the central structure of C(k) since, in the on- 
shell Born approximation, a particle of energy E probes 
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Figure 16. (Color online) Components of the diffusion tensor: 
(solid red line), Db (dashed gray line) and Df (dotted 
blue line) for the 3D incoherent- (top) and coherent- (bottom) 
speckles cases. The dotted lines are power-law fits (Dg oc 
E lu ) to the data in the low and high energy limits. The insets 
show the transport anisotropy factors / D\ and Db/Db • 



C*(k) inside the k-space sphere of radius 2k E centered at 
the origin. In this regime the coherent- and incoherent- 
speckles arc very similar. The most interesting effect 
appears for 2k E a±^ > 3.8 (i.e. E > 1.8E a± ), when the 
bumps of C(k) contribute to scattering and transport. 
The scattering time t e becomes highly anisotropic [see 

Fig. 13(b)] and the orbital dominating is distorted 
compared to the incoherent case. As a result, £)jf is re- 
duced and the corresponding anisotropy factor drops by 
a factor of ~ 4. This effect happens to be strong enough 
to lead to the inversion of the transport anisotropy and 
we find < Dl < D z for E > 1.8E cr± [sec inset of 
Fig. 16(b)]. 



C. Localization 

In order to analyze strong localization effects, we now 
solve the self-consistent equation (77) for the 3D case 
in the long time limit (ui — > 0). A threshold en- 
ergy E c (mobility edge) appears, solution of D^(E C ) = 
det{D B (£ c )} 1/3 = fi/VSnm. For E < E c , one finds 



D*(cj,E) - 0+ - iujtf oc (E) for uj -> 0, where Li oc (E) 
is a real positive definite tensor. As in 2D, it character- 
izes exponential localization within the propagation ker- 
nel (31) with the anisotropic localization tensor h\ oc (E). 
The localization tensor is diagonal in the same basis as 
the Boltzmann diffusion tensor Db. Explicitcly, we have 



(86) 



where Lf^ c = dct{Li oc (£')} 1 / 3 is the unique solution of 




IS 



arctan 



/av 
'B 



(87) 



For E > E c , T)*(lj,E) converges to a real definite pos- 
itive tensor when uj — > 0. It describes anisotropic nor- 
mal diffusive dynamics, characterized by the propaga- 
tion kernel (29) where 'D(E) is replaced by the quantum- 
corrected diffusion tensor 



D»(£) ee lim T)Juj,E) 



l - 



3?rm 2 {Df(E)Y 



D B (£). 



Figure 17 shows the components of Li oc (for E < E c ) 
and D* (for E > E c ) for the single-, incoherent- and 
coherent-speckles cases, and for typical parameters of 
Refs. [57, 58]. As already mentionned in Sec. VIA the 
behavior of Li oc and D* is completely determined by that 
of Db in our approach. The anisotropics of Li\ oc (E) are 
the square roots of those of T)-b(E) [see Eq. (86)] and the 
anisotropics of D*(i?) are the same as those of Db(-E') 
[see Eq. (88)]. Therefore, as for Db, we observe that the 
anisotropy factors of Li oc and D* arc nearly indepen- 
dent of E, except for the inversion of anisotropy of the 
coherent-speckles case. In the single-speckle case we find 
Lf oc /L? f ~ 3.2 and D Z J Dt' y ~ 10. For the incoherent- 
speckles configuration we find L^ z /L^ c ~ 1.3 and 
D'*' z /D^ ~ 1.8. For the coherent-speckles configura- 
tion we find the same values at low energy, and at high 
energy we have Df/D* ~ 0.5 and D z /D* ~ 2.1. Fig- 
ures 17(al), (bl) and (cl) present the results in the local- 
ized regime. At low energy, using the scaling of D^(E) 

obtained previously we find L? oc {E) oc (Dl/Df) 1/2 E 3 / 2 . 
When E increases, L™ ac (E) grows and finally diverges at 
E c . In the diffusive regime [see Fig. f 7(a2), (b2) and (c2)] 
the quantum corrections are significant only close to E c , 
while for higher values of E, D*(E) ~ D B (£'). There- 
fore, in the high E limit we have D${E) cx (D^/D^)E 5 ^ 2 
as found previously (see Sec. VII B). For the coherent- 
speckles case, we recover the inversion of anisotropy pre- 
dicted in Sec. VII B. For the parameters of Fig. 17(c) it 
occurs in the diffusion regime. For higher values of V n , 
however, it can be in the localization regime. 
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Figure 17. (color online) Components of the localization tensor Li oc (left side, for E < E c ) and quantum-corrected diffusion 
tensor D, (right side, for E > E c ) in the single- (upper row; V R = 7.1E a± ), incoherent- (central row; Vr = 0.35E a± ) and 
coherent- (lower row; Vr = 0.35Ecr ± ) speckles cases. The components of Db are plotted for comparison (thin lines on the right 
column). We have used the parameters of Refs. [57, 58] and Fig. 4. 



VIII. ABOUT THE 3D MOBILITY EDGE 

The self-consistent approach used above is expected to 
fairly describe the quantum transport properties [62, 75, 
80]. It gives some quantitative estimates consistent with 
numerical calculations [124] and experimental data [56, 
125]. It however has two main flaws. 

On the one hand, it predicts that, just below the mo- 
bility edge, the localization length diverges as Lf ac {E) oc 
(E c — E)~ u with v = 1 and, just above the mobil- 
ity edge E c , the corrected diffusion tensor increases as 
D%(E) cx (E - E c ) s with s = 1. Those values of the 
critical exponents v and s are consistent with the predic- 
tion s = v[d — 2) of the scaling theory [23, 126] and they 
arc independent of the choice of cut-off that we made. 
However, it is known from advanced numerical calcula- 
tions on the disordered tight-binding model of the An- 



derson model [127, 128] and from experiments [52] that 
they are not correct. The correct value of the critical 
exponents in 3D is v = s = 1.58 ± 0.01 [127, 128]. In 
order to reproduce this value, it seems necessary to take 
into account the fractal nature of the wave functions at 
the critical point [129], which is beyond the scope of the 
self-consistent theory of AL. 

On the other hand, in contrast to critical exponents, 
the location of the critical region (mobility edge, E c ) is 
a non-universal quantity and should be determined from 
microscopic theory. In this respect, the on-shell approx- 
imation is qucstionnablc because it neglects the strong 
modification of the spectral function induced by the dis- 
order. Including this effect is possible within the self- 
consistent theory and we discuss below a new method to 
do it. 
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A. Energy renormalization 

As discussed above, the main failure of the on-shell ap- 
proximation is that it completely neglects the structure 
of the disorder-dependent spectral function, which may 
renormalize energies, and thus strongly affect the value 
of E c . In order to improve the method, one could in prin- 
ciple use the more sophisticated approach of Ref. [124], 
which do incorporate the spectral function, and provides 
values of E c in agreement with numerical calculations in 
the Anderson model. For continuous disorder, one may 
rely on the approach of Rcfs. [60, 64], which has been 
applied to several standard models of disorder. However, 
since we are interested in continuous disordered poten- 
tials with fine anisotropic structures, these methods are 
hardly practicable. From a numerical point of view, es- 
timates of necessary ressources seem out of present-day 
possibilities. In order to overcome this issue, we have 
proposed in Ref. [79] an alternative method based on the 
assumption that the leading term missing in the on-shell 
approximation is the real part of the self energy, 



Z'{E,k) = P 



dk' C(k-k') 



(2n) d E-e* 



(89) 
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Figure 18. (Color online) Comparison of the mobility edge 
as calculated with the SCBA method (the full black squares 
are the results obtained by A. Yedjour and B. van Tiggelen in 
Ref. [64], that we reproduce here), with the on-shell method 
(E' c , red crosses) and with the renormalized self-consistent 
approach (corrected E c , thick blue circles) mobility edges, for 
an isotropic 3D speckle potential. When comparing to Fig. 8 
of Ref. [64], note that in Ref. [64] the reference of energy 
is the minimum value of the disorder and that we have the 
correspondences E^ = E a /2 and U — Fr. 



where P is the Cauchy principal value, see Eq. (43). 
This term renormalizes the energies: A quasi-particle of 
momentum k has an energy E, solution of E — e(k) — 
£'(£■, k) = 0. Here, we incorporate T,'(E,k) into the 
theory sclf-consistcntly and by averaging, in first approx- 
imation, its k-anglc dependence. It amounts to replace 
the on-shell prescription by e(k) = E' = E — A(E) with 



relation function reads [62, 64] 



A(E) = 



1 

4?r 



e(k) = E-A(£) 



dO fc £'(£,k) 



(90) 



Within this approach, all previous quantities are now re- 
garded as functions of E 1 instead of E. It does not change 
the overall energy dependence of the quantities discussed 
above, but may be important for direct comparison to 
energy-resolved experimental measurements. In the fol- 
lowing we concentrate on the 3D mobility edge E c . It 
is the solution of E c — A(E C ) = E' c , where E' c is deter- 
mined using the on-shell approach. The above equation 
is solved self-consistcntly for E c , and A can be regarded 
as an energy shift. 



B. Three-dimensional, isotropic disorder 

In order to validate our approach, we first consider 3D 
correlated disorder with an isotropic correlation function. 
In this case, other methods, such as the self-consistent 
Born approximation (SCBA), can be used to calculate 
the self-energy [60, 64] hence providing a test-bed of our 
approach. For a speckle disorder obtained inside an inte- 
grating sphere lit with a laser beam, the real-space cor- 



C(r) = Vi 



sin (|r|/a) 
(MAx) 2 



(91) 



with a the correlation length. The associated power spec- 
trum (sec appendix B) is isotropic and bears the same 
infrared divergence as the anisotropic models of 3D disor- 
der considered in this work: C(k) cx l/|k| when |k| — > 0. 
It is therefore a relevant model for comparison. Fig- 
ure 18 shows the on-shell mobility edge E' c calculated as 
in Sec. VII C (see also Ref. [62]), the true mobility edge 
E c calculated by our method [self-consistent theory with 
renormalization of the energy, such that E c is the solution 
of E c — £' (E c , He 1 ) = E' c ] , and the mobility edge found 
using the self-consistent Born approximation in Ref. [64] . 
As it is clearly seen in Fig. 18, the disorder-induced modi- 
fication of the spectral function plays a major role for the 
prediction of the mobility edge. While the on-shell mo- 
bility edge, E' c , is positive and increases with energy, the 
corrected mobility edge, E c , as calculated either by the 
method of Ref. [64] or by our self-consistent approach 
with renormalization of energies, is negative and mainly 
decreases with E. In addition, we find that the renormal- 
ized self-consistent approach predicts values of E c in very 
good agreement (within 5 — 7%) with those of Ref. [64]. 



C. Three-dimensional disorder with structured 
correlations 

These results support our method to estimate E c , 
which we now apply to anisotropic disorder in the single- 
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Figure 19. (Color online) Angular dependence of Y*'(E C , fcgyk) 
(thin solid lines) for the single- [upper row; with 6 — (k,k z )] 
and the coherent- [lower row; with O = (k, kz) where the 
different colors represent the variations for k rotating around 
kz] speckles, for different values of Vr (indicated on the fig- 
ure). The horizontal solid blue line is the mean value and the 
dashed blue lines represent the standard deviation around the 
mean, both calculated over the k-solid angle. 
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Figure 20. (color online) On-shell (E' c ) and corrected {E c ) mo- 
bility edges versus the disorder amplitude Vr for the single-, 
incoherent- and coherent-speckles cases and for the parame- 
ters of Fig. 4. 



speckle, incoherent-speckles and coherent-speckles con- 
figurations. The mobility edge is found by searching the 
root of the self-consistent equation (90). Note that the 
averaging of the angular dependence of E' in Eq. (90) 
is justified a posteriori by the weak k-angle variations 
of E' found around its mean value at E c (with standard 
deviations less than 10 — 15%). This is illustrated in 
Fig. 19 which presents the angular variations obtained 
numerically in the calculation of A(E C ) in the single and 
coherent-speckles cases, for typical values of V R . 

The on-shell (E' c ) and renormalized (E c ) mobility 
edges for the three anisotropic models of disorder are 
shown in Fig. 20. As for isotropic disorder, it is eye- 



catching that the shift of the energy states completely 
changes the behavior of the mobility edge. While the 
on-shell mobility edge, E' c , is positive and increases with 
V R , we find that the renormalized mobility edge, E c , is 
negative and decreases with V" R . For V R < E a± , this be- 
havior is qualitatively similar to that obtained in Fig. 18 
(see also Ref. [64]). For larger values of Vr, E c further 
decreases, consistently with the idea that it should ap- 
proach the percolation threshold deep in the classical dis- 
order regime (V R ^> E a± ) [130]. 



IX. CONCLUSIONS 

Disordered potentials with finite-range correlations are 
often characterized by a counter-intuitive and interesting 
behaviour [25, 28-30, 32, 79]. These are directly related 
to the microscopic statistical properties of the potential, 
hallmarked by the disorder correlation function. In this 
paper we have focused on anisotropy effects in 2D and 3D 
correlated disorder. We have quantitatively studied the 
transport and localization of matter waves by using pcr- 
turbativc transport theory [81] and a standard on-shell 
self-consistent approach [75]. The latter, first pioneered 
by Vollhardt and Wolfe [110, 114], remains the most pow- 
erful, quantitative, microscopic approach to Anderson lo- 
calization in dimension higher than one {d ^ 2), in spite 
of the unavoidable problem of describing the physics in- 
side the critical region in d > 2. Within this approach, 
we have characterized incoherent diffusion, quantum cor- 
rected diffusion and localization tensors versus the parti- 
cle energy. We have found rich diffusion and localization 
properties. A striking result is that weak structured cor- 
relations can induce strong anisotropy effects. We have 
supported the general theory with specific examples dis- 
cussing speckle potentials in 2D and 3D. 

2D configuration - In the 2D case, we have considered 
an anisotropic Gaussian correlation function as used in 
Refs. [53, 54]. The energy- dependences of relevant quan- 
tities are studied: For E <C E a± , in the white-noise limit, 
we find t e £ oc 1 for the scattering time and Db oc E for 
the Boltzmann diffusion tensor, which are both isotropic. 
For E > E a± , we find t e £ oc \f~E~ and D B oc E 5/2 . As 
a general rule, the anisotropy of the disorder (£), of the 
scattering time (£ s ) and of Boltzmann diffusion (£b) are 
all different. The scattering time shows an inversion of 
anisotropy from £ s > 1 (for £ > 1) at low energy to 
£ s = l/£ (< 1) at high energy. In contrast, the trans- 
port anisotropy is always £b > 1 (for £ > 1) but shows a 
strongly nonmonotonic behaviour as a function of energy 
with a marked maximum at E ~ E a± . The anisotropy of 
localization is simply the square root of that of transport. 
For typical experimental parameters, we found that it is 
very small in observable regimes, except for very strongly 
anisotropic disorder. So far, experiments have only stud- 
ied the classical regime [53, 54] and our study offers scope 
for future studies of quantum transport and localization 
in 2D speckle potentials. 
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3D configurations - In the 3D case we have discussed 
two configurations recently used to study Anderson lo- 
calization (AL) of matter waves [57, 58]. The energy 
dependence of relevant quantities are the following: For 
E <C E a± , we find t e £ oc 1 and is slightly anisotropic, 
and Db oc E and is significantly anisotropic, which is due 
to anisotropic suppression of the white-noise limit in the 
models we used. For E E a± , we find t e £ oc \f~E~ 

and Db oc E 5 / 2 , both being anisotropic. Again, the 
anisotropics of the different quantities arc all different. 
In the single- speckle case, we found almost constant 
anisotropy factors {D z jD x ' y ~ 10 and L\ JL X ^ C ~ 3.2), 
and experimental data can be compared to these pre- 
dictions almost independently of the energy distribution. 
Although no precise value has been extracted from the 
experiment of Rcf. [57], the data indicate significantly 
larger anisotropy. Further analysis would be required to 
clarify the origin of such discrepancy. In the incoherent- 
speckles case we also found almost constant anisotropy 
factors (D^ Z /DX ~ 1.8 and L z oc / L* o v c ~ 1.3). In the 
coherent- speckles case, bumps in the disorder power spec- 
trum induce a strong inversion of anisotropy as a function 
of energy. The anisotropy factor measured in Ref. [58] 
was shown to be in fair agreement with the theory. The 
inversion of the transport anisotropy was however not ob- 
served because the images were taken in the (y, z) plane. 
It only gave access to D« = D Y and D z = [D x + D z )/2, 
which do not show the inversion. In order to observe it, it 
is required to image the atoms along the transport eige- 
naxes and to tunc the balance between the populations 
of low- and high-encrgy states. 

Mobility edge in 3D - We have further studied the 
behaviour of the 3D mobility edge. To do so, we have 
extended the on-shell approach and proposed a way to 
renormalize energies. We have found a striking agree- 
ment of our method with the more involved method 
based on SCBA developed in Ref [64] for isotropic dis- 
order. The effect of renormalizing energies does not al- 
ter the overall energy dependence of te quantities dis- 
cussed above, but may be important for direct com- 
parison to energy-resolved experimental measurements. 
It is worth discussing our predictions in view of what 
has been experimentally achieved so far. Comparing 
to Ref. [57] (single- speckle configuration), our calcula- 
tions significantly differ from experimental values (e.g. for 
V R = 600 nK x k B ~ 7.1E a± , we find E c ~ -300 nK x k B 
while +900 nK x k B is measured). However, the method 
used in Ref. [57] to infer E c from the localized frac- 
tion neglects the disorder-induced distortion of the en- 
ergy distribution. It is questionable because the latter is, 
in particular, necessary to account for negative energy 
states (i.e. below the disorder mean value). Comparing 
to Rcf. [58] (coherent- speckles configuration), we find that 
A(£' c ) as calculated here is of the same order of magni- 
tude as the heuristic shift introduced in Rcf. [58] (e.g. for 
V B = hx680Uz~0.35E a± , we find A{E c )/h = -390 Hz 
and the heuristic shift is —225 Hz). A precise test of 
the present theory would however require a reliable de- 



termination of the energy distribution in ultracold-atom 
experiments, which is not available so far. 

Finally, our results and method may provide a guide 
line to future experiments investigating the so-far un- 
explored effect of anisotropy in quantum transport of 
matter waves. In the case of ultracold atoms, to which 
our study directly applies, the transport properties can 
be probed by direct imaging of the atoms and control 
of the energy. The effects discussed in this manuscript 
can be expected for other kinds of waves and/or other 
models of disorder, and arc particularly relevant to new 
systems where the disorder correlations can be con- 
trolled [10, 57, 58, 67, 131-133]. 
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Appendix A: Intensity kernel 

In this section we show the step-by-step calculation 
of the long-time and large-distance limit of the intensity 
kernel given by Eqs. (56), (57) and (59) and the diffusion 
tensor Eq. (61). 

As explained in Sec. II D, the solution of the Bethe- 
Salpetcr equation (20)-(21) can be obtained by inverting 
the operator A = l-G<g)Gt U [see Eq. (23)]. To this aim, 
we diagonalize the operator G®G^ U in the (q, a;) = (0, 0) 
limit. We thus solve 

f dk' 

J W fEM ' = A ^> k (A1) 

where U^ k , = t/ k ,k'(q = 0,a; = 0,E) and f E ^ = 
G(£,k)Gt(£,k) [see Eq. (25) for q = and uj = 0]. 

1. Preliminary remark 

First, let us notice that we have 

fEM = ^^-A(E,-k), (A2) 

where A(E, k) is the spectral function defined in Eq. (12) 
and r s (i?,k) is the scattering mean free time defined in 
Eq. (14). 
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2. Properties of Eq. (Al) 

The main properties of Eq. (Al) and of its eigcnfuc- 
tions are listed below: 

1 . The eigenvalues X E and the eigenvectors (f> E k of 
Eq. (Al) are real. 

Proof. By multiplying Eq. (Al) by Gt(.E, k), we obtain 



Hk' 

^ M k B k , Gt (E, k') <f>% M , = X E Gt (E, k) d>% M , 



where -Aijf k , 



(A3) 

G*{E, k) E/£ k , G(£, k'). The latter is 

Hermitian since G^{E,V)* = G(E,k.) and U£ w is real 
and symmetric. Therefore all the eigenvalues X E are real. 
By taking the complex conjugate of Eq. (A3), dividing 
by G(E, k) and comparing it to Eq. (Al), we obtain that 
the functions <fi E k are real. 

If U^ k , is positive-definite, the eigenvalues X E are 
positive. In particular, this is always true in the 
Born approximation [134]. When U£ k , is symmet- 
ric and positive-definite, we can write it as U^ k , = 

I (i^ ( 5k,k"rfk"Q k '' ;k ', where d k » > and Q is an or- 
thogonal operator. For any vector of components Xk, 
we have / £ )d ffj^ M^av = / ^ k |y k | 2 > 0, 
where y k = J -{0jd G^(E, k')x k 'G/ k < k . It shows that 
Af k 3 k , is positive definite. Its eigenvalues X E are therefore 



positive. 



□ 



2. The eigenvectors <f> E k can be chosen to satisfy the 
orthonormalization condition 



(A4) 



Proof. This is an immediate consequence of the fact that, 
according to Eq. (A3), the functions G^(E, k) ^ k are 
eigenfunctions of the Hermitian operator M k 5 k , . □ 

3. The eigenvectors <f> E k satisfy the completeness re- 
lation 

fEM £ <t> n E,K = W 5 ( k ~ k ')- ( A5 ) 



Proof. This follows from the fact that the eigenfuntions 
G$(E, k) <j> Et . of the matrix Mjf k ,, Eq. (A3), form a com- 
plete basis. □ 



Proof. We multiply both terms of Eq. (Al) by <j> E k , and 
sum over n. Equation (A6) is recovered by using the 
completeness relation Eq. (A5). □ 

5. The most important property of Eq. (Al) is that 
one of the eigenvalues is 

AF 1 = 1, (A7) 

and the corresponding eigenvector is proportional 
to the inverse scattering mean free time: 



^Tk 1 = ^ 



'e(^k)]- 



%A(E,k) [T.(£,k)]-i 



(A8) 



(27T) 



Proof. This is a direct consequence of the Ward identity 
[110]: 



AE k (q,w,£) = 



dk' 

(2n) a 



U k! u>(q,u),E) AG k (q, W ,£), 

(A9) 



where AE k (q,w,£) = E(B + ,k+) - St(B_,k_) and 

AG k (q,w,£0 = G(£+,k + )-Gt(£_,k_). For (q,w) = 
(0, 0) it becomes 

/dk' 
p^v/^AE k (0,0 1 E). (A10) 

When comparing Eq. (A10) to Eq. (Al), we obtain that 
AE k (0,0,£;) = -ih/r s (E,k) is a solution of Eq. (Al) 
with unit eigenvalue. Using Eq. (A2) and the or- 
thonormalization condition (A4) one then easily finds 
Eq. (A8). □ 

6. The eigenfunctions <j) E k have the parity properties: 



/,n=l _ ±n=\ 



IJ E,-k 



U EM 



for n > 1. 



(All) 
(A12) 



Proof. This is a consequence of the parity of the vertex 
^kk" in particular, U\ _ k , = Ug v . Using Eq. (A6) 
we have £„ ^e <^,k = E„ A^.-k <^_ k , , which 
can only be satisfied if the eigenfunctions <fi E k have a 
well defined parity. The eigncfunction 4> E = £ is given by 
Eq. (A8) and it is even. In addition, using Eqs. (A2) and 
(A8) in the orthonormalization condition (A4), we have 
J A(E,k)(j) EM = for n > 1. Which shows that 



(27T)<" 



are odd functions of k. 



□ 



Solution of the BSE 



4. The irreducible vertex function U^ k , can be ex- 
pressed as 



^ B k' = E A 



n aji in 
E fEM 9 EM! 



(A6) 



Note first that, if Eq. (Al) could be diagonalizcd with 
all eigenvalues different from one {X E ^ 1 for all n), it is 
straightforward to show, using Eq. (A5), that we would 
have A^(0,0,E) = £ n [l/(l - \%)]fk<%<%, . In this 
case no diffusion would be observed. As noticed above, 
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however, the conservation of particle number, through 
the Ward identity, imposes that there is one eigenvalue 
equal to one. As there is no other conserved quantity in 
the system we are considering, we can assume that the 
eigenvalue A = 1 is not degenerated and that there is 
a finite gap between this eigenvalue and the rest of the 
spectrum when (q, uj) — , [135, 136]. This suggests the 
following ansatz for the solution of the BSE (20)-(21) [see 
Eq. (23)], in the small (but non-zero) q and uj limit: 



Then, making use of the parity properties of the func- 
tions <^ fc [Eqs. (All) and (A12)], t e £ (even func- 
tion of k) and V k e(k) (odd function of k), we fi- 



nally obtain ^(q, u, E)f E .M = ^7k{q, E)/ x /ti(r \ 



<J> fn „ P1 -f ^(q, W ,£)^(q,w,£) 



A(q,w,£7) 



1 



^ 1-A 1 



EM' 



(A13) 

where </> k (q, w > -^0 an< ^ 1 + -Ml; w 7 are solutions of the 
eigencquation 

/Air' 
J^~ d U^f^^^^E) 

= [l + X(fl,u,E)]4i(q,u,E). (A14) 

The latter are the first eigenvalue and eigenvector at 
small (q, uj), and reduce to Eqs. (A7) and (A8) when 
(q, uj) = (0, 0), respectively. We then write /k(q, uj, E) = 
Jem + -Fic(q, uj, E) the expansion of / k (q, w, E). Making 
the ansatz <p^(q,uj,E) = V a„(q, uj, E)cj> E k , we find 



\(q,uj,E) = J2 



a n {q,uj,E) f dk 



ai(q,uj,E) J (2tt) 

(A15) 

Finally, the coefficients a„(q, ui,E) are found by im- 
posing that Eq. (A13) solves the BSE. After some al- 
gebra one finds ai(q, uj,E) = 1 and a„(q, uj,E) = 



l-A' 



dk 

(27T)<> 



b° Ek F k (q,uj,E)cf, Ek , forn>l. 



4. On-shell approximation 



We now proceed to the on-shell (weak disorder) ap- 
proximation, and we neglect the effect of disorder on the 
spectral function. Equation (A2) becomes 



f E M^^A (k,E), 



(A16) 



where t e k is the on-shell scattering mean free time 
[r £k = r a (E,k E k)}, A (k,E) = 2tt5[E - e(k)] and e(k) 
are, respectively, the disorder-free particle spectral func- 
tion and dispersion relation. An explicit calculation of 
the small (q, uj) expansion of fk(q,uj,E), gives [137] 

( ' 2 



F k (q,L0,E) 



EM 
H 2 



[fed - q • V k e(k)] 



2r 3 - r 3 - 

+ huj [q • V k e(k)] - -§£ [q. V k e(k)] 2 



ft 3 

xA {k,E) + O(uj 2 ,q 3 ,q 2 uj). (A17) 



^4 jk F k (q,cd,£) ^,k- 



B,k' 

where 7 k is given by Eq. (60) and A(q, uj,E) = 
2N a (E) [itvui - Hq-B(E)-q] /Hr~\) with the diffusion 

tensor of Eq. (61). The solution of the BSE is thus 
given by Eq. (56) with Eqs. (57) and (59). Note that 
this expression for the diffusion constant is quite general 
(only the on-shell approximation has been made), pro- 
vided that the full irreducible vertex function U is con- 
sidered in the eigenequation (Al). In Sec. VA the Born 
and Boltzmann approximations are made U = Ub [see 
Eq. (62)]. 



Appendix B: Isotropic disorder 

For disorder with isotropic correlation function, we 
define, as in Ref. [62], p(k,9) = C{k\k - k'|) = 
C(2k\ sin(6*/2)|) , where 8 is the angle between the unit 

vectors k and k and k = |k| = |k'|. In this case, rotation 
invariance ensures that the eigenproblem (62) is solved 
by cylindrical (2D) or spherical (3D) harmonics. 

1. Two-dimensional case 

In the 2D isotropic case, inserting the cylindrical har- 
monics Zq = 1, = cos(W) and Z^ 1 = sin(Z0) into 
Eq. (62), we find 



J^d6p(k E ,9) cos(W) 



\ 1.171 



d0p(k E ,i 



(Bl) 



where I > and m € { — 1, +1} are integer numbers. In 
particular, we find X l E ° = 1 in agreement with Eq. (A7). 
They are doubly-degenerated for I > and the corre- 
sponding normalized eigenfunctions are proportional to 
the orthonormal cylindrical harmonics, with the prefac- 
tor determined by the normalization condition (A4): 



and 



b l=0 - 



h E.k 



Zo(d) 




z±\e)\ 



d6>p(k E ,0' 



^ d0' p(k E ,6>) 



(B2) 



(B3) 



In the calculation of the diffusion constant, it is actually 
possible to see that only the first term plus the I = 1 
terms (with m = — 1, +1) in the summation of the right- 
hand side of Eq. (61), contribute to the diffusion coeffi- 
cient. More precisely the on-shell scattering mean free 
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time t e £ does not depend on k, v x (respectively v y ) is 
a 27r-periodic and even (resp. odd) function of 9, and 
(resp. is 27r/Z-periodic and even (resp. odd). 

Therefore, when performing the angular averaging of the 
product t e itVi&g ^ i n Eq. (61), one finds that only the 
term with I = 1 and m = +1 (resp. m = —1) couples to 
v x (resp. v y ) and contribute to (resp. D B ). Then, 
inserting Eqs. (Bl), (B2) and (B3) into Eq. (61), we find 



D B (E) = 



HE 



mN (E) J 2 * d9 (1- cos 6) p(k E , 6) 



(B4) 



This formula agress with the result of Ref. [62] , obtained 
by a different approach. 



E=6.3xl0 3 E„ E=6.3xl0 'e„ E=63E n 
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2. Three-dimensional case 

In the 3D isotropic case, proceeding in a similar way, 
we find that the eigenvalues of Eq. (62) are given by 



f*dO sin 9 p(k B , 9) P t (cos( 
f*d0 &mOp(k E ,0) 



(B5) 



with the index I = 0,1, +oo and m = —I, —I + 1, +1 
and where Pi (cos 6) are the Legendre polynomials. The 
eigenvalues are (21 + l)-degenerated and the correspond- 
ing normalized eigenfunctions are propotional to or- 
thonormal spherical harmonics, with the prefactor de- 
termined by the normalization condition (A4): 



em 



Y, r 




2tt J d6' sm9'p(k B ,9'), (B6) 



In the calculation of the diffusion constant, using the 
same type of symmetry arguments as in the 2D case, we 
find that only the I = 1 (with m = — 1, 0, 1) terms couple 
to v and contribute in the summation of Eq. (61). We 
thus find 



D B (E) 



hE 



1 



3tt mN (E) /* dO sin (9 (l - cos (9) p(k B ,9) ' 

(B7) 

which agrees with the expression found in Ref. [62]. 



3. Three-dimensional isotropic speckle 

A simple model of 3D speckle with isotropic correlation 
properties, is found when considering the light pattern 
obtained inside an integrating sphere lit by a laser beam 
of wavevector k L . The real-space correlation function is 
given in Eq. (91) and the associated power spectrum 



V 2 n 2 cr 2 
C(k) = »„ , 0(2tr- 



(B8) 



is isotropic. Although this isotropic model is unrealistic 
from an experimental point of view, it is useful here in 



Figure 21. (Color online) Isotropic 3D speckle. Eigenvalues 
of Eq. (62) (top row) for the isotropic 3D speckle with power 
spectrum given by Eq. (B8). Topology of the main eigenvec- 
tors 4>% k contributing to D§ (bottom row), D| (2 nd row) and 

-Dg (3 rd row) [with the parametrization k = (ka;,ky,k z ) = 
(sin 9 cos <f>, sin 6 sin <f>, cos 9)], the red lines locate the nodal 
lines. From left to right E = 6.3 x lCr :i £ CT , E = 6.3 x KT 1 ^ 
and E = 63E a . 



two respects. First, it bears the same divergence as the 
anisotropic 3D models of disorder considered in Sec. Ill: 
C(k) oc l/|k| when |k| — > 0. Second, several proper- 
ties of this model are analytical and known [61, 62], and 
therefore provides a test our numerical methods. 

As done previously, for the diagonalization of the inte- 
gral operator (62) we use 2 7 x 2 7 points regularly spaced 
on the k-space shell |k| = k E . Some eigenfunctions and 
eigenvalues of Eq. (62) are presented in Fig. 21. We in- 
deed find spherical harmonics [see Eq. (B6)], and the 
eigenvalues A^, agree well with theory [Eq. (B5) with 
C given by Eq. (B8), not shown on the figure]. We 
further incoporate these results in Eq. (61). Figure 22 
presents the numerical results for the Boltzmann diffu- 
sion constant (red dots) which agree very well with the 
analytic formula (solid black line) found when incorpo- 
rating Eq. (B8) into Eq. (B7). Note that we recover the 
same asymptotic behaviours as for our anisotropic cases: 
D B (E) oc E for E/E a < 1/2 and D B (E) oc E 5 ' 2 for 
E/E a > 1/2. In particular, those tests show that the 
discretization used here correctly treats the |k| — > di- 
vergence. 



Appendix C: Conductivity 

1. Einstein relation 

As presented in Sec. HE, we expect <t(uj = 0) oc D in 
the linear response regime. Here we calculate <r B (w = 0) 
in the Boltzmann approximation and verify this relation 
explicitly, which enables us to find the proportionality 
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Figure 22. (Color online) Boltzmann diffusion coefficient 
for the isotropic 3D speckle configuration of power spectrum 
given by Eq. (B8). The solid black line is the theoretical 
prediction, red dots are numerical results. 



factor in Eq. (34). 

Let us first rewrite the Boltzmann diffusion tensor as 



where Jk is the rcnormalizcd current vertex 
J k 2tt ^ A 



T = v + T ^ T 



(C2) 



We want to calculate the conductivity <x B in the ladder 
approximation. We have to evaluate 









— ~ 


H* (C3) 

M — 1 









where T is defined in Eq. (55). It reads 
dk 



(27T) 

dk dk' 



£ VifE,k v j 



(C4) 



(27r) d (2tt) 



3 w»/B,krk 1 k'(o,o, J5)/B )k /^-. 



As 1^,(0,0,25) = Ex^iirSf^t^ [138]. and 
/s,k — t s £.-<4.o(22, k)/7i, one eas ily finds 



i / n 27T f / \ 27T 



(C5) 



El 



X'' 



Therefore, we have cr B = 27riVo(-E , )DB/ / S. We have thus 
verified Einstein's relation for the classical dc conductiv- 
ity in anistropic disorder. 

2. Current vertex renormalization 

The DC conductivity <r B in the Boltzmann approxima- 
tion reads (see appendix CI) 



(C6) 



where Jk, the renormalized vertex function, is given by 
Eq. (C2). Diagrammatically we can absorb this renor- 
malization in one of the vertices as shown in Eq. (C7). 
This is a standard procedure for anisotropic scattering, 
which is presented for example in Ref. [6]. 



(C7) 



3. Weak-localization correction 

a. The cooperon 

We calculate the bare cooperon correction, with renormalized current vertices, Diag. (68) translates into 



(C8) 



Considering that the dominant contribution in the integral comes from Q ~ k + k' ~ [see Eq (73)], and that 
Je k ~ 2(r B ^/h) 3 Ao(E, k) in the on-shell approximation [139], we get 



iai'L(u,E) = 2 —— ( 

<- x r ' 7 hN (E) \ 



dQ 



ft 2 E ' k / J (27r) d + HQ ■ T>b(E) ■ Q 



(C9) 
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b. Hikami contributions 



We now calculate the Hikami corrections [see Diags. (69) and (70)] 

(CIO) 

In the same way as before, and using the on-shell approximation formulas [139] G(E, k)/^;^ ~ ~^( t e k/h) 2 Ao(E, k) 
and G*(E,k)f E ,k ~ i{r Eic /h) 2 A (E,k), we get A<x (ffl) ~ &(T(h 2 ) and 

A^)(^£) = Aa^ i) (u>,E)+Aali 2) {io,E) = (Cll) 

_2 ,Jkji /• dk' Jk^> /-_dQ 1 

HNq(E) h Te *J (2ir) d BkM ' lEM ' h 'J (2ir) d —iHu) + HQ ■ D B (E) ■ Q ' 



c. Corrected conductivity tensor 

We now consider the quantity Jk — / j^ppc t^Bk,k'/ B k'^k'- Using the relation C/ B k.k' = ^e^e k^s k'' an< ^ 

the parities of the functions 0™ ^ [see Eqs. (All) and (A12)], one can show that 

dk' 

f/B k ,k'/s,k'Jk' = hv. (C12) 



(27r) d 

Therefore the Hikami contributions renormalize one of the Jk/h back to the bare vertex v, and we have 

2 ^ dQ ^ (C13) 



V ' 7 W V ' ; W V ' ' HNo(Ey h 3 B ' k/ 7 (27r) d -i/ia; + ftQ ■ D B (£) • Q 

which gives the final expression (74). 
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